ODC  FILE  COPY,  AD  A 0 4 2 2 6 3 


AFGL-TR-77-0044 


MRDA 

A MEDIUM  RESOLUTION  DATA  ANALYSIS  CODE 
FOR  THE  HP2100  COMPUTER 


D.  Kryger  and  D.  Robertson 


AERODYNE  RESEARCH,  INC. 
Bedford  Research  Park 
Crosby  Drive 
Bedford,  MA  01730 


C8/^r  AVA!i;.8LE 
PERMIT  FL'Li  y LI 


TO  COG  e;>s  i:ot 

'’.•pK  CJ'i 


January  1977 


Final  Report  for  Period  26  March  1976  to  31  Decenober  1976 


Approved  for  pidsllc  release;  distribution  unlimited 

The  views  and  conclusions  contained  In  this  document  are  those  of  the  authors 
and  should  not  be  Interpreted  as  necessarily  representing  the  official  policies, 
either  expressed  or  Implied,  of  the  Defense  Advanced  Research  Projects  Agency  or 
the  U.S.  Government. 

This  research  was  supported  by  the  Defense  Advanced  Research  Projects  i\gmicy 
of  the  Department  of  Defense  iinder  ARPA  Order  2656,  monitored  by  U.S.  Army 
Missile  Command  under  MIPR  A31687-63-7257. 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
HANSCOM  ABF,  MASSACHUSETTS  01731 


D D C 

»j?ar?nn/7r?i 


ETETEinnSl 

E 


ARPA  Order  No.  2656 

Program  Code  No.  6E20 

Contractor:  Aerodyne  Research,  Inc. 

Effective  Date  of  Contract 
26  March  76 


Contract  No.  F19628-76-C-0173 

Principal  Investigator  and  Phone  No. 

Dr.  James  Draper  617-275-9400 

AFGL  Project  Scientist  and  Phone  No. 
Dr.  Brian  P.  Sandford  617-861-4281 

Contract  Expiration  Date:  31  Dec.  76 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
Documentation  Center.  All  others  should  apply  to  the  National 
Technical  Information  Service. 


Unclassitied 


SeCURlTV  CLASSIFICATION  OF  Th*^  *»A0‘  O*  » 

/ j)  REPORT  DOCUMENTATION  PAGE 


/* 


( 


K BfcPB"'''  »ft<unEn  

' AFGL-TR-77-(i<^44  U' 


\2  COVT  ACC  essics  NO 


4.  Title  (m%d  Submit) 


MRDA  - A Medium  Resolution  Data  Analysis  Code  / 
P"or  the  HP  2100  Computer  , — __5T  . 


7.  AUTMORr»> 

D./Kryger 
D./ Robertson 


RF/ID  INSTRUCTIONS 
HKhCKi:  COMPl.ETINO  FORM 


3 PFCl'^  Ff  T'S  catalog  NUMBER 


S.  T\rt  Or  Ff  PONT  4 PERIOD  COVERED 

Final  Hoport 
;i/26/7G  - 12/31/76 

Sr  PERE  QNMlN^i  Qftfi  REPORT  NUMBER 


Ari-RR-97 


9 PER^OPMING  organization  NAME  AND  ADDRESS 

Aerodyne  Research,  Inc. 

Crosby  Drive,  Bedford  Research  Park 
Bedford,  Massachusetts  017.'10 


B CONTRACT  OrVranT  NUMBERf») 

i 1;U9()2877(>-C-017:1/>-^'-— 

' ' -A*'/?  F /r  , - , — L ' 

Iment.  p'rojI'IIT.'  T«I! 


AI  EA  A ^'R<  UNIT  numbers 

^ — r /f  — 

62101 F/ 76  7|li  101  -7  JJ. 

ARPA  ORDER  2656  / 


11,  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Geophysics  Laboratory 
Hanscom  AFB,  Massachusetts  01731 
Monitor  - B.  Sandford/OPR 


'Ll 


REPORT  DATE 

Ja 


77  / 


13  NUMBER  n<=  PAGES 

85 


<A.  MONITORING  AGENCT  NAME  » AOORESSfll  dlllerml  from  Controlfmit  OllUfJ 

U.S.  Army  Missile  Research  & Development  Command 
DARPA  Projects  Office  (DRDMI-NS) 


IS  security  class,  fof  thfi  r..port; 

Unclassified 


Redstone  Arsenal,  AL  .35809 


16.  DISTRIBUTION  STATEMENT  fof  Ihim  W»POfll 


Ml?'  nip 

7^..; 


A- 


Approved  fdt  public  release;  distribution  unlimited 


W.  DISTRIBUTION  STATEMENT  (ol  tht  sbMtrmct  tntertd  in  Block  20,  It  dlUtrtnt  from  Report) 


18.  supplementary  NOTES 


This  research  was  supported  by  The  Defense  Advanced  Research  Projects  AgengV' 
of  the  Department  of  Defense  under  ARPA  Order  No,  2656,  monitored  by  Army 

r- 


Missile  Command  under  MIPR  A31687-63-7257 


19  KEY  WORDS  (Conttmjo  on  roeorte  eido  It  nmeeatmry  mnd  Idonltty  by  bfocfc  number} 

Atmospheric  Transmittance 
Minicomputer 


I ABSTRACT  (Continue  on  revere#  tide  If  ntcettory  mnd  Identity  by  ti)rk  numbtr) 
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I 


INTRODUCTION 


The  task  of  generating  the  Medium  Resolution  Data  Analysis  (MRDA)  code  was  under- 
taken to  provide  medium  resolution  predictions  of  atmospheric  transmission  and  radiation. 
Equally  important  is  the  requirement  that  the  software  operate  on  the  HP2100  minicomputer. 
This  enables  a direct  convenient  inclusion  of  atmospheric  effects  in  the  radiation  studies  being 
carried  out  at  AFGL.  MRDA  is  basically  a combination  of  the  high  resolution  code  HITRAN^^^ 
and  the  low  resolution  code  LOWTRANS.^^^  The  basic  input -output  structure  of  MRDA  is 
similar  to  that  of  LOWTRAN3,  with  additions  to  improve  its  spectral  resolution. 

The  MRDA  code  is  now  operational  on  the  HP2100  mini -computer.  Because  of  severe 
limitations  on  core  size  the  code  is  divided  into  seven  segments.  The  MRDA  code  correctly 
reproduces  the  spectral  structure  of  atmospheric  transmittance  calculations  when  compared 
to  results  using  HITRAN.  However,  MRDA  tends  to  overestimate  the  absorption  at  hi^er 
altitudes  in  the  spectral  region  near  the  center  of  a stron^y  absorbing  line.  The  framework 
for  calculating  atmospheric  radiation  is  built  into  the  program,  althou^  these  calculations 
are  not  done  in  the  present  version  of  MRDA. 

Section  II  outlines  the  instruction  format  for  using  MRDA.  In  Section  III  a discussion  of 
the  organization  of  MRDA  is  presented  along  with  examples  of  how  various  types  of  calculations 
are  performed.  Section  IV  discusses  the  capabilities  of  the  code  and  Its  accuracy.  A listing 
of  the  computer  code  is  given  in  Appendix  A,  supplemented  by  a flow  chart  for  each  segment 
(Appendix  B)  and  a definition  of  symbols  (Appendix  C).  A listing  of  the  program  used  to  gen- 
erate the  absorption  coefficients  is  in  Appendix  D. 


(1)  R.  A.  McClatchey  et  al. , "AFCRL  Atmospheric  Absorption  Line  Parameters," 
AFCRL-TR-73-0096,  January  1973,  /C 

(2)  J.  E.  A.  Sebly  and  R.  A.  McClatchey,  "Atmosi^ieric  Transmittance  from  0.25 
to  28.5  pm:  Computer  Code  LOWTRAN3,"  AFCRL-TR-0255,  May  1975. 


1 1 INSTRUCTION  FOR  USING  MRDA 


Since  MRDA  is  very  close  to  LOWTRAN3  in  concept  and  operation,  tlie  injjuts  to  the  two 
profframs  are  virtually  identical.  For  this  reason  much  of  the  description  presented  here  is 
a paraphrase  of  the  LOWTRAN3  operation  instructions.  The  input  medium  for  MRDA,  instead 
of  a card  deck,  is  a disk  JOB  file.  This  file  contains  all  the  input  information  that  would 
normally  appear  in  the  input  card  deck  for  LOWTRAN3  with  a few  changes.  A listing  rf  the 
JOB  file  (called  MDATA)  is  shown  in  Figure  1.  The  many  different  atmospheric  data  contained 
in  LOWTRAN3  are  not  used  by  MRDA.  MRDA  uses  none  of  the  spectral  atmospheric  data  and, 
in  order  to  conserve  memory,  uses  only  one  model  atmosphere.  Presently,  the  file  MDATA 
contains  the  1962  U.S.  Standard  Model  Atmosphere. 

In  general,  it  is  only  necessary  to  change  the  last  four  lines  of  the  JOB  file  (referred 
to  here  as  input  lines  1-4)  in  order  to  run  the  program  for  a given  problem.  The  formats 
for  the  last  four  input  lines  and  their  application  will  be  discussed  next.  The  user  has  the 
option  of  different  paths  (horizontal,  slant,  vertical,  to  space)  and  the  option  of  reading  in 
radiosonde  data  to  replace  the  model  atmosphere. 

A.  INPUT  DATA  AND  FORMATS 


The  data  necessary  to  specify  a given  problem  are  given  on  the  last  four  lines  of 


the  JOB  file. 

They  are  as  follows: 

Input  Line  1 

MODEL,  IHAZE,  ITYPE,  LEN,  JP,  IM,  NLDAT,  IRAD 

FORMAT  (713) 

Input  Line  2 

HI,  H2,  ANGLE,  RANGE,  BETA,  VIS 

FORMAT  (6F  10.3) 

Input  lane  3 

VI,  V2,  DV 

FORMAT  (3F  10.3) 

Input  Line  4 

DVM 

FORMAT  (F  10.2) 

Definitions  of  the  above  quantities  will  be  discussed  in  Section  B. 

If  the  quantity  Mt)DEL  given  on  input  line  1 is  set  equal  to  0 or  2 (which  means 
meteorologlcai  data  arc  used  as  in[)Ut  to  the  program),  then  the  above  card  sequence  (and 
format  for  input  line  2)  is  changed.  These  cases  will  be  described  in  Section  C. 
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Figure  1.  (cont.)  Kxample  of  the  MRDA  JOB  File 


B.  liASIC  INS'nUU'TIONS 

The  various.  QUantitioe  to  ho  spi'clfied  on  oaclt  of  Uic  four  control  oanla  are 
discussed  in  tins  section 

1.  Input  Line  1 MODEL.  IHAZE,  ITYPE,  LKN.  .IP,  IM,  NLDAT,  '.HAD 

The  parameter  MODEL  either  selects  the  model  atmosphere  or  specifies  that 
meteorological  data  are  to  be  used  in  |)lace  of  the  standard  model.  IHAZE  sjxjcifies  whether 
aerosol  attenuation  is  to  be  included  in  the  calculation  or  not.  For  any  problem  the  atmos- 
pheric path  must  lie  specified  as  one  of  three  types  according  to  ITYPE  and  LEN.  The  rest 
of  the  quantities  given  on  input  line  1 (which  can  be  left  blank  if  not  required)  provide  the 
user  with  options  to  suppress  printing  (JP),  and  to  input  a new  model  atmosphere  (IM, 
NLDAT).  The  options  for  the  above  parameters  and  their  uses  are  stated  and  described  in 
detail  below; 

MODEL  = 0 if  meteorological  data  are  specified  (for  horizontal  paths  only)* 

= 1 selects  1962  US  STANDARD  Model  Atmosphere 
= 2 if  a new  model  atmosphere  (or  radiosonde  data)  is  to  be  instered.* 

MRDA  currently  uses  only  the  1962  U.S.  Standard  Model  atmosphere.  If  any  of  the  other 
MIWTKAN.J  motlel  atmospheres  are  desired  a new  JOB  file  can  be  created  with  the  ap- 
pn>pnate  model  atmosi)here  in  place  of  the  U.S.  Standard  and  setting  MODEL  = 1. 

MRDA  will  still  identify  the  model  as  the  1962  U.S.  Standard.  Since  MRDA  assumes  the 
earth's  radius  is  6371.23  km  this  may  introduce  minor  differences  when  a model  other 
than  the  U.S.  Standard  atmosphere  is  used, 

IMAZE  = 0 means  no  aerosol  attenuation  included  in  tlie  calculations, 

IHAZE  = 1 or  2 if  aerosol  attenuation  is  required  (see  also  Input  Line  2). 

If  IHAZE  is  set  equal  to  1 or  2 and  visual  range  (VIS)  is  not  specified  on  Input  Line  2,  then  the 
program  will  automatically  select  visual  ranges  of  23  km  or  5 km,  respectively. 


♦ In  these  cases  the  format  for  Input  Line  2 changes  (see  nonstandard  conditions)  Section  C 


ITYPK  = 1 for  a horizontal  (constant  pressure)  patli. 

= 2 for  a vertical  or  slant  path  between  two  altitudes. 

= 3 for  a vertical  or  slant  patli  to  s))ace. 

The  TYPE  1 path  should  not  be  confused  with  a lon5>,  path  where  the  local  height  at  tlie 
end  of  the  trajectory  is  slKnificantly  different  from  the  beginning  heij^t.  In  such  a case, 
spxjcify  the  path  according  to  rrYPE  = 2 (see  Section  BZ) 

LEN  = 0 for  normal  operation  of  program. 

LEN  = 1 selects  the  downward  TYPE  2 patli  shown  in  Hgure  2(e). 

The  parameter  LEN  can  be  ignored  (that  is,  left  blank)  for  the  majority  of  cases.  It  need 
only  be  used  for  a downward  looking  path  (H2  > HI)  when  two  paths  are  possible  for  the 
same  input  parameters.  In  such  a case,  a computer  printout  statement  will  be  given 
indicating  that  the  user  has  two  choices  for  the  problem  and  that  the  shorter  path  (see 
Figure  2<d))  has  been  executed.  Set  LEN  = 1 for  the  longer  case . 

JP  = 0 for  normal  operation  of  program. 

JP  = 1 to  suppress  printing  of  transmittance  table 
IM  » 1 when  radiosonde  data  are  to  be  read  in  initially 

IM  = 0 for  normal  operation  of  program  or  when  subsequent  calculations  are 
to  be  run  with  MODEL  = 2 

NLDAT  = Number  of  levels  to  be  read  in  for  MODEL  = 2 

Note  that  IM  and  NLDAT  are  only  used  when  MODEL  = 2. 

IRAD  is  a flag  to  indicate  whether  or  not  a radiation  calculation  is  to  be  performed 
and  is  defined  as  indicated  below. 

IRAD  = 1 perform  radiation  calculation. 

IRAD  = 0 do  not  perform  radiation  calculation. 

Only  the  IRAD  * 0 option  should  be  used  in  the  present  version  of  MRDAl  In  the  case 
of  IRAD  = I an  additional  input  is  required  between  the  first  and  second  input  lines.  This 
input  is  as  follows: 

EMBS,  TRACK  FORMAT  (2F  10.3) 


(d)  (•) 


Figure  2.  Geometrical  I’atli  Conliguration  for:  (a)  Horizontal  Paths  (Type  1),  (b)  Slant 
Patlis  Between  Two  Altitudes  Ml  and  H2  (Type  2),  and  (c)  Slant  Paths  to  Space  (Type  3). 
For  downward  looking  paths  where  HMIN  < H2  < HI,  two  trajectories  are  possible  as 
indicated  in  (d)  and  (e).  From  Reference  (2). 
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EMISS  and  TBACK  define  the  emissivity  and  temperature  of  a background  radiation 
source  located  at  the  beginning  of  the  atmospheric  path.  TBACK  is  in  units  of  °K. 

In  the  case  where  MODEL  = 2,  the  new  atmosphere  (model  or  radiosonde  data)  Is 
inserted  between  input  Lines  1 and  2 (see  Section  C), 

2,  Input  Line  2 ill,  H2.  ANGLE,  RANGE,  BETA,  VIS 

Input  Line  2 is  used  to  define  the  geometrical  path  parameters  for  a given  problem. 

HI  = Initial  altitude  (km) 

H2  = Final  altitude  (km) 

ANGLE  = Initial  zenith  angle  (degrees)  as  measured  from  HI 

RANGE  = Path  length  (km) 

BETA  = Earth  centre  angle  subtended  by  HI  and  H2  (degrees) 

VIS  =Sealevel  visual  range  (km) 

Figure  2 shows  the  different  atmospheric  paths  that  are  possible.  It  is  not  neces- 
sary to  specify  every  quantity  given  above;  only  those  required  to  adequately  describe  the 
problem  according  to  the  parameter  ITYPE  (as  described  below)  are  needed. 

(1)  Horizontal  Paths  (ITYPE  = 1) 

(a)  Specify  HI,  RANGE  and  VIS  only 

(b)  If  nonstandard  meterorological  data  are  to  be  used,  that  is,  if 

MODEL  = 0 on  Input  Line  1,  then  the  following  parameters  must  be  specified  on  Input 
Line  2:  H1,P,  T,  RH,  WH,  WO,  VB,  RANGE  according  to  FORMAT  (3F  10.3, 

2F  5.  1,  2E  10.3,  2F  10.3),  where  P,  T,  DP,  RH,  WH,  and  WO  are  the  pressure  (mb), 
temperature  (°C),  dew  point  temperature  (°C),  relative  humidity  (%),  HgO  density  (gm  m ^), 
respectively.  It  is  necessary  to  specify  all  of  the  quantities  underlined  with  a full  line 
and  one  of  the  quantities  underlined  with  a dashed  line. 

(2)  Slant  Paths  to  Space  (ITYPE  = 3) 

specify  Hi,  ANGLE  and  VIS 

(3)  Slant  Paths  Between  Two  Altitudes  (ITYPE  = 2) 
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(a)  SiKJoify  HI,  H2,  ANill.l!,  imd  VIS;  or 

(b)  Siwilfy  Ml,  AN(11,I';,  IIANCK,  iind  VIS;  or 

(c)  S|)ecify  HI,  H2,  RANr.E,  and  VIS. 

For  cases  (b)  and  (c),  the  program  will  calculate  112  and  ANGLI*]  respectively, 
assuming  no  refraction  and  then  proceed  as  for  case  (a).  This  method  of  defining  the  pro- 
blem should  be  used  when  refraction  effects  are  not  important,  e.g. , for  ranges  of  a few 
tens  of  km  at  zenith  angles  less  than  80°,  It  can  also  be  used  for  larger  angles  (including 
90°)  provided  that  the  path  lies  within  one  atmosjjheric  layer.  Slant  paths  to  space 
(ITYPE  =3)  or  between  two  altitudes  (ITYPE  =2)  in  wdiich  both  ANGLE  and  RANGE  are 
not  defined,  are  not  possible  in  MRDA.  In  LOWTRAN3  these  paths  require  the  calculation 
of  the  initial  zenith  angle  of  the  path  by  use  of  subroutine  ANGLE.  MRDA  does  not  contain 
this  routine: thus  these  paths  are  not  allowed.  If  this  case  is  attempted,  an  error  message 
will  be  printed  and  execution  halted. 

In  the  case  where  MODEL  =2  the  new  model  atmosphere  (or  radiosonde  data)  is 
inserted  between  Input  Lines  1 and  2 (see  Section  c). 

3.  Input  Line  3 VI,  V2,  DV 

The  spectral  range  over  which  transmittance  data  are  required  and  the  spectral 
increments  at  which  the  continuum  calculations  are  to  be  made  are  given  by; 

VI  = initial  frequency  in  wave  numbers  (cm 

V2  = final  frequency  in  wave  numbers  (cm  where  V2  > VI 

DV  = frequency  increment  (or  step  size)  (cm  ^)  for  continuum  calculations. 

4 -1 

(Note  that  v = 10  /X  where  v is  the  frequency  in  cm  and  X is  the  wavelength  in  microns, 
and  that  DV  can  only  take  on  values  which  are  a multiple  of  5cm  \) 

4.  Input  Line  4 DVM 

The  variable  DVM  defines  tlie  spectral  steps  at  which  the  medium  resolution  calcul- 
ations are  done  and  the  wavenumber  interval  of  tlie  final  output.  A value  of  0.05  cm~^  should 
be  used,  because  it  agrees  with  the  sjxjctral  resolution  of  the  library  tape. 


c. 


ATMOSPHKIUC  MODEL 


rhnn.'  options  arc  available  if  atmospheric  transmittance  calculations  are  required 
for  nonsUindard  conditions.  Here  nonstandard  refers  to  conditions  other  than  those  specified 
by  the  model  atmosphere  in  MDATA.  The  three  options  enable  tlic  reader  to  insert: 

(H  His  own  model  atmosphere(s)  in  place  of  tlie  standard  model,  provided  that 

the  data  are  in  exactly  the  same  format  and  are  si)ecified  at  the  same  altitudes 
as  the  latter.  In  this  case  the  appropriate  print  statements  in  MRDA  (that 
identify  the  atmospheric  model  used)  must  be  changed  correspondintdy. 

(2)  An  additional  atmospheric  model  (MODEl;=2),  which  can  be  in  the  form  of 
radiosonde  data.  The  data  need  not  be  specified  at  the  same  altitudes  as  in 
the  standard  models. 

(3)  Meteorological  conditions  for  a given  horizontal  path  calculation  (MODEL  = 0 case). 
This  was  discussed  in  Section  B-2  and  will  not  be  repeated  here. 

The  first  of  these  options  requires  the  most  effort  and  needs  no  further  discussion 

here. 

1.  Additional  Atmospheric  Model  (MODEL  =2) 

As  stated  in  Section  B-2,  a new  model  atmosphere  can  be  inserted  between  Input 
Line  1 ard  2 provided  the  parameters  MODEL  and  IM  are  set  equal  to  2 and  1 respectively 
on  Input  Line  1.  The  number  of  atmospheric  levels  to  be  inserted  (NLDAT)  must  also  be 
specified  on  Input  lane  1.  The  appropriate  meteorological  parameters  and  format  for  the 
atmospheric  data  are  given  below. 

Z,  P,  T,  DP,  RH,  ^i,  WO,  AHAZE  (FORMAT  (3F  10.3,  2F  5.1,  2E  10.3,  2F  10.3)) 


where 

Z 

= altitude  (km) 

P 

= pressure 

T 

= ambient  temperature  (°C) 

DP 

= dew  ix)int  temperature  (°C) 

RH 

= relative  humidity  (%) 

WH 

= water  vapor  density  (gm  m 

WO 

= ozone  density  (gm  m ) 

AHAZE 

= aerosol  number  density  (cm 

■i 


Note  that  it  is  only  necessary  to  specify  those  quantities  underlined  with  a full  line  and 
either  of  the  quantities  underlined  witli  the  dashed  line. 

If  the  aerosol  nu_rnber  density  was  not  measured  as  a function  of  altitude  and  tlie  reader 
wishes  to  include  aerosol  attenuation  in  the  calciJation,  set  IHAZK  = 1 on  Input  line  1,  In 
this  case  MRDA  will  use  the  aerosol  models  already  contained  in  tlie  program  and  interjjolate 
to  give  aerosol  number  density  values  at  the  same  altitudes  as  the  radiosonde  (or  new  model 
atmosphere)  data.  The  program  \«11  then  look  for  a sea  level  visual  range  (V15)  to  be 
specified  on  Input  Line  2.  If  VIS  is  not  specified,  a 23  km  sea  level  visual  range  will  be 
assumed.  If  aerosol  attenuation  is  not  required,  set  IHAZE  = 0 on  Input  Line  1 as  before. 

D,  RUNNING  MRDA 

As  already  mentioned  all  the  inputs  for  MRDA  are  contained  within  the  JOB  file 
MDATA.  Before  running  MRDA  the  HP2100  system  should  be  in  the  normal  keyboard  mode 
since  the  JOB  file  uses  FMGR.  To  execute  MRDA  simply  type 

♦RUN,  JOB,  MD,  AT,  A 

Before  actually  running  MRDA  tw'o  things  should  be  done.  First  mount  the  appropriate 
MRDA  input  tape  containing  the  spectral  extinction  coefficients  corresponding  to  the 
frequencies  of  interest  on  UNIT  0,  Second,  depending  on  the  value  of  IRAD,  MRDA  uses 
1 or  2 disk  files  (TRANl  and  TRAN2).  If  the  previous  MRDA  run  ended  other  than 
normally,  one  or  both  of  these  files  may  still  exist  in  the  directory  and  should  be  purged 
before  ■running  MRDA  once  again.  The  output  from  MRDA  will  appear  on  the  line  printer. 


♦Any  JOB  file  with  the  structure  of  MDATA  (possibly)  containing  a different  model  atmosphere 
could  just  as  wll  be  used. 
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III  MRDA  DESCRIPTION 

A.  DESCRIPTION  OF  MRDA  SOFTWARE 

MRDA  is  desiKned  to  make  medium  resolution  atmospheric  transmission  cal- 
culations in  the  range  of  1800  to  6000  cm  over  a wide  variety  of  geometrical  paths.  These 
calculations  are  carried  out  within  MRDA  in  three  stages.  The  first  part  of  the  calculation 
consists  of  predicting  the  continuum  transmissions  (H  O,  N , molecular  scattering)  along  with 
aerosol  absorption  for  the  chosen  path  and  wavenuinoer  interv'al.  These  calculations  are 
carried  out  by  what  is  essentially  LOWTRAN3  with  all  the  spectral  calculations  removed. 

In  the  process  of  computing  the  continuum  results,  intermediate  values  are  saved  for  later 
use.  For  example  the  pressure,  temperature  and  altitude  of  each  layer  in  the  geometric 
path  are  stored.  In  addition,  if  a radiation  calculation  is  requested  (IRAD  = 1)  the  trans- 
mission throit-^  each  layer  is  also  stored  in  a disk  file  named  TRANl.  Finally  the  atmos- 
pheric concentration  of  H.,0  and  0„  in  each  layer  along  with  the  molecular  density  (of  all 
gases)  for  the  particular  path  through  each  layer  is  computed  and  saved. 

In  the  second  part  of  MRDA  the  medium  resolution  spectral  calculations  are  performed. 
The  atmospheric  information  needed  to  carry  out  these  calculations  have  been  transfered  from 
part  1.  'rhree  magnetic  tapes  contain  a complete  library  of  extinction  coefficients  necessary 
to  compute  the  medium  resolution  results.  These  tapes  are  generated  using  the  hi^ 
resolution  transmission  program  (HlTRAN).  The  tapes  are  organized  in  10  cm  ^ blocks 
over  the  1800  to  6000  cm  ^ range  that  MRDA  operates.  In  each  wavenumber  block 
exctinction  coefficients  for  6 molecular  species  at  9 pressure-temperature  points  are  defined 
at  wavenumbers  which  depend  on  the  structure  of  the  absorption  spectra  for  the  particular 
species  in  that  wavenumber  block.  For  each  species  the  wavenumbers  were  chosen  so  as  to 
define  the  spectra  by  identifying  the  strongest  lines  in  the  region.  The  extinction  coefficients 
are  then  calculated  at  the  line  center,  points  slightly  removed  from  the  line  center  and  mid- 
way between  adjacent  strong  lines. 

Using  the  information  from  part  1,  the  library  tape  data,  and  possibly  the  disk  file 
TRANl,  the  program  calculates  the  total  spectral  transmission  for  the  geometric  path 
and  frequencies  defined  by  tJie  input.  In  order  to  obtain  the  extinction  coefficient  from  the 
tape  at  the  particular  pressure,  temperature  and  frequency  required,  the  program  per- 
forms linear  interpolations  over  the  pressure -temperature  matrix  defined  on  the  tape  and 
then  over  frequency.  Once  having  calculated  the  total  transmission  at  a particular 
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frequency,  that  frequency,  the  transmission  and  radiation  results  are  written  to  a disk  file 
TRAN2. 

Part  3 of  the  MRDA  software  combines  tlie  continuum  and  medium  resolution  results 
The  continuum  results  from  part  1 are  linearly  interpolated  to  obtain  transmist.ion  values  at 
the  medium  resolution  frequencies.  Finally  the  total  transmission  is  computed  at  the 
medium  resolution  frequencies  and  written  to  a disk  file  for  later  use. 

B.  PROGRAM  ORGANIZATION 

Structurally  MRDA  is  composed  of  7 separate  software  segments  or  overlays. 

As  described  in  Section  A the  7 segments  can  be  combined  into  3 well  defined  parts. 

Part  1 (modified  LOWTRAN3)  consists  of  segments  1 through  5.  Part  2 is  segment  6 and 
part  3 segment  7.  A complete  listing  of  all  segments  and  subroutines  is  shown  in 
Appendix  A.  Appendix  B contains  a segment  by  segment  detailed  logic  flow  diagram 
for  MRDA.  Appendix  C has  a description  of  the  variables  used  in  each  segment. 

To  improve  the  readability  and  ease  in  following  the  logic  with  MRDA,  a particular 
notational  convention  for  indices  has  been  followed  throughout.  This  convention  is  outlined 
below: 

Index  variables 

1 - frequency  index 

L - atmospheric  level  number.  0 to  NL 

LL  - layer  number,  numbered  sequencially  along  the  geometric  path 
K - species  index 

In  addition  many  variables  beginning  with  1,  L,  or  K refer  to  quantities  related  to  these 
indicies . 


t I 
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C.  MRDA  I.LHUAHY  TAPES 


The  MHDA  I.ibrary  tapew  conUiin  tin-  speclnil  ab8orf)tlon  coefficients  for  the  six 
atmospheric  species  which  have  significant  absorrrtion  in  the  1800  - 0000  cm”'  region. 

The  species  are;  H20»  CO.,,  O^,  N^O,  CO  and  Cli^.  The  absorption  coefficients  are 
calculated  at  the  selected  (P,  T,  u)  points  and  then  written  onto  a tape  that  is  accessed 
by  MRDA,  The  CDC6600  computer  at  AFGL  was  used  to  generate  the  MRDA  tapes. 

1.  Choice  of  Spectral  Absorption  Coefficients 

Because  of  the  severe  constraint  on  the  amount  of  available  storage  on  the 
HP2100,  the  data  in  the  MRDA  library  tape  are  organized  so  as  to  define  the  absorption 
spectra  for  tlie  species  in  as  compact  a form  as  possible.  Thus  nine  pressure-tempera- 
ture (P,  T)  points  are  used  to  describe  the  atmosphere,  and  the  total  number  of  wave- 
numbers  points  within  each  block  is  limited  to  250. 

The  choice  of  the  (P,T)  points  is  based  on  the  expected  range  of  atmospheric  (P,  T) 
values.  Figure  3 shows  the  (P,  T)  profile  of  the  U.S.  Standard  Model  Atmosphere,  the 
estimated  limits  of  (P,  T)  variability,  and  illustrative  radiosonde  data  taken  from  several 

3 

Af'CL  Teal  Ruby  Missions.  The  heavy  dots  within  the  circles  show  the  nine  (P,  T)  points 
at  which  the  spectral  absorption  coefficients  are  calculated.  Pressure-temperature  points 
for  pressures  below  100  mb  are  not  included  in  the  tape  at  this  time  due  to  the  problems 
with  the  spectral  structure  of  the  transmittance  at  these  hippier  altitudes. 

A considerable  savings  in  the  total  number  of  wavenumber  points  at  which  the  spec- 
tral absorption  coefficients  must  be  stored  is  obtained  by  Identifying  the  stronger  spectra] 
lines  within  each  wavenumber  block.  When  one  or  more  of  the  species  have  no  strong 
lines  within  a block,  the  absorption  coefficients  are  calculated  at  only  four  points.  After 
locating  the  strong  peaks,  the  extinction  coefficients  are  calculated  at  the  peak,  at  points 
±0.05  cm  removed  from  the  peak,  and  at  a point  half  way  between  two  adjacent  peaks. 

The  spectra  from  one  strong  line  to  the  next  are  described  by  five  points.  Between  these 
points  the  sf)ectral  absorption  coefficients  are  obtained  by  linear  interpolation. 


3.  B.Sandford,  et  al.,  "Aircraft  Signatures  in  the  Infrared  1.2  to  5.5  Micron  Region," 
AFGL-TR-76-0133,  Air  Force  Geophysics  Laboratory  (OPR),  Hanscom  AFB,  Mass. 
01731,  (Juno,  1976) 
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>Fit?uro3  Temperature  and  pressure  variations  of  the  atmosphere. 

Radiosonde  data  from  several  Teal  Ruby  Missions  ai*e  indicated. 
I The  outer  lines  indicate  the  approximate  range  of  atmospheric 

r temperature  fluctuations.  The  center  line  is  the  U.S.  Standaixi 
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2 .  Calculation  ^ the  Spectral  Absorption  Coefficienta 

The  spectral  absorption  coefficients  are  calculated  from  the  molecular  line 
parameters  on  tlie  AFCL  IHTHAN  Data  Tape.  ^ A listing  of  the  proKram  is  ^iven  in 
Appendix  D.  The  data  for  the  input  parameters  are  given  by  the  following  sequence  of 
five  cards: 

1.  NI>TPTS  (E) 

2.  P (1).  I = 1,  NPTPTS  (8K10.0) 

3.  T (I),  1 = 1,  NPTPTS  (8E10.0) 

4.  W (M),  M = 1,  7 (7E10.3) 

5.  VI,  VV2,  DV,  BOUND  (6F10.3) 

The  input  quantities  are : 

NPTPTS  = number  of  (P,  T)  points 
P = pressure  values 

T = temperature  values 

2 

W = species  column  density  = 0.2H9E20  molecules/cm  for  (li„0,  CO^, 

O3,  N2O,  CO,  CH^,02> 

VI  = lower  frequency  limit,  cm 

VV2  = upper  frequency  limit,  cm  ^ 

DV  = frequency  increment,  fixed  at  0.  05  cm 

BOUND  = frequency  from  line  center  beyond  which  a line  is  not  included,  cm 

A value  of  20  cm  ^ is  used  for  the  parameter  BOUND.  Since  oxygen  does  not  have  any 
important  absorption  bands  below  6000  cm  , it  is  not  included  as  one  of  the  six  species 
in  MRDA. 

A Lorentz  line  shape  Is  used  to  calculate  the  contribution  of  each  line  to  the  spec- 
tral absorption  coefficient.  The  absorption  coefficient  at  each  value  Is  a weighted  ave- 
rage ef  three  mono<AroiiMitic  calculations  within  a spectral  interval  of  width  0,06  ozn 
centered  at  v^. 
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IV  DISCUSSION 


The  atmospheric  transmission  code  MUDA  is  o}»rational  on  the  HP2100  minicomputer. 
The  software  has  been  written  in  a format  which  allows  the  user  considerable  flexibility  in 
his  choice  of  atmospheric  conditions  and  paths  for  transmittance  calculations.  In  the  regions 
with  absorbing  molecular  lines,  MRDA  reproduces  the  correct  spectral  strucutre  of  the  atmos- 
pheric transmittance  and  generally  gives  good  quantitative  predictions.  However,  as  discussed 
below,  values  of  the  transmittance  calculated  by  MRDA  for  atmospheric  paths  where  the 
pressure  is  less  than  600  mb  is  low  in  the  immediate  region  of  a strong  absorbing  line.  This 
indicates  that  further  study  into  improving  the  parameterization  of  the  extinction  coefficients 
at  higher  altitudes  is  needed.  The  structure  of  the  HP2100  software  code  is  designed  with 
sufficient  generality  so  that  changes  in  the  parameterization  of  the  extinction  coefficients  are 
quite  easily  incorporated. 

A,  LOW  RESOHmON  TRANSMITTANCE 

The  low  resolution  transmittance  is  calculated  in  the  same  manner  as  LOWTRAN3. 
Included  are  the  contributions  due  to; 

• 4pm  HgO  continuum, 

• 5pm  N.,  continuum, 

• Aerosol  Scattering,  and 

• Molecular  Scattering 

Transmittance  calculations  by  the  MRDA  code  have  been  compared  to  those  calculated  by 
LOWTRAN3  and  have  been  found  to  be  in  exact  agreement. 


B.  ATMOSPHERIC  PATHS 

MRDA  is  designed  to  handle  a variety  of  atmospheric  paths; 

• Paths  to  space,  horizontal  paths,  and  slant  paths  (upward  or  downward). 

The  identification  and  treatment  of  these  paths  are  similar  to  those  in  LOWTRAN3.  However, 
the  high  spectral  resolution  part  of  the  calculation  requires  that  the  absorber  amount  for 
each  atmospheric  layer  be  considered,  not  just  the  total  equivalent  sea  level  absorber 
for  the  entire  path.  In  MRDA,  these  absorber  amounts  arc  calculated  in  segments  three. 


A: 
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four  and  five  and  stored  in  the  matrix  WW.  They  are  passed  on  to  segment  six  in  the 
matrix  W(3AS,  Because  the  concentrations  of  water  vapor  and  ozone  vary  independently 
of  the  concentrations  for  the  well  mixed  gases,  they  are  also  calculated  and  passed  on  in 
the  matrices  W1120and  W03,  respectively.  These  quantities  are  calculated  and  stored  for 
all  atmospheric  paths. 

C.  HIGH  UESOLUTION  TRANSMITTANCE 

The  transmittance  due  to  the  spectx*al  structure  pf  absorption  by  the  atmospheric 
gases  is  calculated  in  segment  six.  Basically,  the  calculation  depends  on  two  different 
inputs:  the  amount  of  an  absorbing  species  along  a segment  of  the  path  and  the  extinction 
coefficient  for  that  species  at  the  temperature  and  pressure  of  the  segment.  The  total 
transmittance  is  then  the  product  of  the  transmittances  for  all  the  segments.  The  results  at 
each  wavenumber  are  both  printed  out  and  stored  on  disk  for  later  use. 

For  verification  of  the  high  resolution  transmittance  calculations,  comparisons  were 

made  to  calculations  made  using  the  HITRAN  data  tapes.  The  HITRAN  calculations  were 

done  every  0.01  cm~^  and  degraded  to  a resolution  of  0. 1 cm  ^ by  using  a triangular  slit 

function  of  half-width  0.  1 cm  Horizontal  paths  were  used  that  had  a constant  amount  of 
20 

absorbing  species  (10  molecules).  The  transmittance  was  calculated  at  several  different 
temperatures  and  pressures  in  the  spectral  interval  2300-2310  cm  , which  has  strong  CO2 
absorption. 

At  lower  altitudes  the  agreement  between  transmittances  calculated  by  HITRAN  smd 
.MRDA  was  good.  The  regions  of  both  strong  and  weak  absorption  were  calculated  properly 
by  MRDA.  The  agreement  was  generally  to  within  10%;  when  the  transmittance  is  very  small 
(less  than  0. 10),  the  difference  would  become  much  larger.  This  is  because  the  line  shape 
in  not  sampled  as  extensively  in  MRDA  as  in  HITRAN.  As  seen  in  Figure  4 MRDA 
reprcKfuoes  the  overall  spectral  structure  of  the  transmittance  calculated  by  HITRAN. 

At  higher  altitude  the  agreement  between  MRDA  and  HITRAN  is  not  as  good  because 
MRDA  overestimates  the  absorption  due  to  strong  lines.  Figure  5 shows  the  ttr:msmittance 
calculated  by  MRDA  and  HITRAN  for  a horizontal  path  at  500  mb  pressure  and  2G0°K  temper- 
ature. Figure  6 show.s  a similar  transmittance  at  100  mb  and  220°K.  This  indicates  that 

further  investigation  into  the  calculation  and  storage  of  the  spectral  absorption  coefficients 
is  required.  The  initial  choice  of  the  extinction  coefficients  was  determined  by  the  requirement 


Transmittance 


P"i(4ure  f>.  Comparison  of  MRDA  and  HITRAN  Transmittance  Calculations 
at  500  mb  Pressure 
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Wavenumber,  cm 

Figure  6.  Comparison  of  MRDA  and  HTl’RAN  Transmittance  Calculations 
at  100  mb  Pressure 
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that  the  number  of  entries  on  the  data  tapes  not  be  too  large  and  that,  at  the  same  time,  the 
correct  spectral  structure  be  reproduced.  The  spectral  structure  of  the  transmittance 
is  calculated  correctly  at  all  pressures  and  temperatures  but  the  absorption  is  overestimated 
(by  one  or  two  orders  of  magnitude  in  some  cases)  when  the  full  widths  of  stron4y  absorbing 
lines  become  much  narrower  than  the  intended  resolution  of  MRDA,  0.  1 cm"^.  Further 
investigation  is  necessary  to  determine  an  approximation  scheme  which  will  quantitatively 
parameterize  these  narrow  lines  better  and  stay  within  the  size  restrictions  of  the  HP2100. 

D,  DOPPLER  BROADENING 

The  spectral  structure  of  the  transmittance  at  altitudes  above  18  km  depends  on 
Doppler  broadening  of  the  absori)tion  lines.  The  Lorentz  width  of  a spectral  line  decreases 
linearly  with  the  pressure  as  one  goes  to  higher  altitudes  and  lower  pressures.  In  this 
framework  Doppler  broadening  can  be  thoug^it  of  as  providing  a lower  limit  to  tJie  linewidth. 
Doppler  broadening  is  approximated  near  line  centers  by  making  the  extinction  coefficient 
independent  of  pressure  for  altitudes  above  18  km.  Since  the  gas  densities  are  very  low  at 
these  high  altitudes,  this  approximation  should  not  significantly  affect  the  total  transmittance. 
However,  quantitative  verification  must  wait  for  a more  accurate  parameterization  of  the 
absorption  coefficients  on  the  MRDA  tape. 

E. ‘  RECOMMENDATIONS 

The  MRDA  software  code  is  capable  of  calculating  the  transmittance  for  a range  of 
atmospheric  paths.  The  code  is  operational  on  the  HP2100  computer.  Because  of  the 

HP2100's  severe  limitation  on  core  size,  it  consists  of  seven  overlays.  As  the  quantitative 
accuracy  of  the  code  needs  improvement,  it  is  strongly  recommended  that  further  inves- 
tigation be  conducted  to  determine  a better  parameterization  for  the  absorption  coefficients. 
The  difficulty  in  getting  the  MRDA  code  operational  on  the  HP2100  has  necessitated  that 
the  major  contractual  effort  be  directed  towards  writing  the  code.  The  result  is  that  further 
effort  is  required  to  improve  the  code's  quantitative  predictions. 

Some  specific  recommendations  for  upgrading  the  code  and  making  it  more  flexible 
are  listed  below. 

• Absorption  Coefficients:  Improve  the  parameterization  of  the  absorption 
coefficients  for  pressures  below  600  mb 


Raciiation:  Incorporato  Uu.-  calculaliun  of  atirjosjiheric  radiation.  Most 
of  the  framework  and  the  necessary  parameters  are  already  incorjxj rated 
in  the  MUD  A code 

Model  Atmospheres;  Increase  the  user's  choice  of  model  atmospheres 
from  the  LJS  Stiuidard  to  include  the  otlier  five  model  atmospheres  in 
LOWTUAN3. 

Input;  Allow  path  description  to  be  made  through  the  teletype  as  well  as 
input  files. 
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COMMON  MM  - MOLEL - 1 1 - Tz - TM I N - P 1 - P2 - PM  I N 

INTEGER  HZ^2-2:.' 

DRTR  NRME  -2HSE  - 2HGT,  .3H3 

DRTR  HZa  - 1 2H23  - HZ''  1 1 2).  2HKM-' - HZ 2 - 1 ,'/2H  5 ■ - HZ ( 2 - 2 > ■ "SHL  M 

PPOGRRM  MR HR  CRLCOlRTES  THE  TPRNSM I TTRNCE  OF  THE  RTMOSPHERE 
FROM  1800  CM-1  TO  6000  CM-1  ':  1 . 54  TO  5.56  MICRONS)  RT  0.O5  CM-l 
3PECTRRL  INTERVRLS  ON  R LINERR  NRVENUMBER  SCREE. 

CRLCULRTION  DEFINED  B'f'  FOUR  CRRD  SEQUENCE  RT  END  OF  JOB  FILE 

C RRD  1 MODEL • I HRZE - I T VPE - LEN » JP  5 I M - NLDRT 
CRPD  2 H1)H2-RNGLE-RRNGE-BETR,VIS 
CRRD  3 VIjVZ-DV 
CRRD  4 DVM 

M0DEL=1  SELECTS  U.S.  STRNDRRD  MODEL  RTMOSRHERE 

MODEL=0  FOR  HORIZ.  PRTH  NHEN  METEOROL.  DRTR  USED  INSTERD  OF  CRRD  2 
PERD  HI  - P(MB  ' - T' DEG  C.'^DEN  RT . TEMP ':  DEG  C'jXREL  HUMIDITY- H20  DENS!'^ 
GM.CM-3:'-  03  DENS  I TY  L GM . M-3 ) » V I S ':  KM ) , RRNGE ':  KM  > NITH  FORMRT  429. 
M0DEL  = 2 NHEN  NEN  MODEL  RTMOSRHERE ':  E . G . RhDIOSONDE  DRTR  J USED. 

DRTR  CRRDS  RRE  RERD  IN  BETNEEN  CRRDS  1 RND  2-  RND  SHOULD  CONTRIN 
RL  T I LUDE  ( KM  - PRESSURE  - TEMR  - DEN  PT  . TEMP  - REL . HUM  I D I TY  - H20  DENS  I TY  - 
0 3 DENSITY. REROSOL  NO.  DENS  I TY ( CM-3 ) RCCORDING  TO  FORMRT  429. 

NOTE  THRT  EITHER  DEN  PT.  TEMP.  OR  REL.  HUMIDIT'i'  CRN  BF  USED 

IF  IHRZE=0  NO  REPOSOL  SCRTTEPING  IS  COMPUTED 

IHhZE=1  IF  REPOSOL  RTTENURTION  REQUIRED  'THIS  IS  USED  IN 

CONJUNCTION  NITH  VI  SURE  PRNGECSEE  CRPD  2 ■ 

IHRZE=1  OR  2 RLSO  GIVE  REPOSOL  RTTENURTION  FOR  23  KM  RND  5 KM  VIS. 
HRZE  MODELS  RESPECTIVELY  IF  VIS=0  ON  CRRD  2 

ITYPE=1.2  OR  3 INDICRTES  THE  TYPE  OF  RTMOSPHERIC  PRTH 
1T','PE  = 3.  VERTICRL  OR  SLRNT  PRTH  TO  SPRCE 
IT  ,'PE  = 2.  VERTICRL  OR  SLRNT  PRTH  BETNEEN  ! NO  RLTITUDES 
ITYPE=1.  CORRESPONDS  TO  R HORIZONTRL  - CONSTRNT  PRESSURE:.'  PRTH 

H 1 =OBSER VER  RL T IT UDE  ■, KM ) 

H2=S0URCE  RLTITUDE  (KM  ' 

RNGLE=ZENITH  RNGLE  RT  HI  (DEGREES) 

RRNGE=PRTH  LENGTH  'IKM) 

BETR=ERRTH  CENTRE  RNGLE 
VIS=VISURL  RRNGE  RT  SER  LEVEL  (KM) 
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IF  ITYPE=1  REflIi  HI  HNI)  RRNGE  IF  ITVPE=:i  READ  HI  RHU  ANGLE. 

IF  ITYPE=2  READ  HI  AND  TNG  OTHER  PARAMETERS  E.G.  H2  AND  ANGLE) 

V1--=1NITIRL  FREQUENCY  •.WAVENUMBER  CM-1  • INTEGER  VALUE 
V2  = FINRL  FREOUENCY  aJRVENUMBER  CM-D  INTEGER  VALUE 
ri'-,'  = FREQUENCY  INTERVALS  AT  WHICH  TRANSMITTANCE  IS  PRINTED 
NOTE:  DV  MUSI  BE  A MULTIPLE  OF  5 CM-1 


DVM  = MErHUM  RESOLUTION  F REQUENC'i'  INTERVAL  'USE  0.05  CM-1) 

IF  'IRAD.NE.O'  I USER  I R CARD  FOR  EMISSIVITY  AND  TEMPERATURE 
OF  A BACKGROUND  BLACF BODY  SOURCE-  INSERT  AFTER  CARD  1. 

EM  I SS= ATMOSPHERIC  EMISSITIVITY  FOR  RADIATION  CALCULATIONS 
TBACK=BACKGROIJND  ?r<te^T*,USED  FOR  RADIATION  CALCULATIONS 
IPAD  =0  ATMOSPHEPnrTADIATION  NOT  CALCULATED 
=1  ATMOSPHERIC  PADIATION  CALCULATED 
ONLY  RIJN  PRESENT  VERSION  WITH  IRAD  =0  ** 
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GO  T 0 1 2 ? 2 - 20O  - 3O0 I PRM 
12  IPRM  = 1 
BET  = 0.0 
BETA  = 0.0 
IXY  = 0 

READ  <5-400'  lATM.NL 

READ  <5- 401)  < HZ  1 ' L ) - L==  1 - 34 ) 

READ  ''5-401^  ':H22‘ L>-L=l-5) 

DO  1 1=1 -HL 

1 READ  <5-402)  2<I  '-P<I  .'-TU  '-WH' : P-WO'::i  ' 

R'E AD  'v  5 - 43 1 ) < V X ' I ' f C7 ' I ) - C7 A < I ) - I = 1 - 44  ' 

READ  <5-405)  <C4‘  I .:• » 1 = 1 - 133) 

READ<5-404)  (CSC  I -- 1 = 1 - 15) 

PI=2.0*flSIN<l .0) 

CA=PI.^180. 

IP=0 

2 CONTINUE 
RE=6371.23 

C JP  NE  0 SUPRESS  PRINT 

READ  <5.400  :'  MODEL.  IHA2E-  ITYPE- LEN- JP  - IM-NLDAT-  IRAD 
WRITE  <6  >400-'  MODEL  - I HAZE  - I TYPE- LEN-  JP  - IM-NLDAT  - IRAD 
200  M=MODEL 

IF'  IRAD. EQ. 0 ' GO  TO  360 
READ<5-406)  EMISS-TBACK 
WRITE <6 >351)  EMISS-TBACK. 

351  FORMAT ' ' .•  " BACKGROUND  SOURCE : EMI SS I V I T Y = " - F5 . 2 > 

15X, '•TEMPERATURE="F5.  1-"  DEGREES  KELVIN"  .•' 

360  IF'- IXY.GT.3)  GO  TO  8 

I F ■ M . EC! . 2 . AND . I M . NE . O ) GO  TO  4 
IF'M.EQ.0)  GO  TO  4 

300  READ<5-406)  HI > H2- ANGLE  - RANGE- BETA- VIS 

WRITE <6- 425)  HI  - H2- ANGLE- RANGE- BETA- VIS 
X1=RE+H1 

IF'  ITYPE.  EU.S.-  GO  TO  560 
IF ' ITYPE. EQ. 1 - GO  TO  8 
X2'=REfH2 

IF' RANGE. EQ.0)  GO  TO  5 

WPITE<6-428)  HI -HZ'- ANGLE -RANGE -BETA- VIS 
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01  52 

0133 

0134 

0135 

3 

0136 

0137 

0138 

4 

0139 

0140 

0141 

0142 

0143 

0 144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0 1 58 

0159 

0160 

500 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

520 

0 1 69 

0170 

0171 

0172 

540 

0173 

0174 

0175 

c 

0176 

c 

0177 

0178 

560 

0179 

0180 

0181 

580 

0182 

0183 

5 

0184 

0185 

c 

0186 

c 

0187 

c 

SEGFl  F IN4  CUMPILER;  HP24177 


I ^ hi' . F.i ' . y . hN [t . hUGL  E . ME. . U ■'  ijO  TO  3 

hMGLE -RCOS^  y . • Hi-Hl  1 . ■t-Xl  "Xl  >.-PhMGE  KhNGE  XI  •>  Ch 
Oh  TO  I 

Xw  = SuP  T ''v X 1 PhMGE  + PfiMGL-'  X 1 + 2 . 0*C0S  '■  RMGL  E"  >'  F<  ' ♦K 1 ♦PhNGF  ■' 
H1=X2-PE 
GO  TO  7 
C DMT  I HUE 

IF'.MLDHT.LE.n  'NLlihT  = ] 

DO  540  L=1.MIDHT 
MHhZE' L '=0. 0 

I F ' N . EC! . 0 ' REHIi ' 5 < 419  ' H 1 ■ P ' 1 > , IMP  • DP  ^ F'H  - NH ' L ' » WO  C L « V 1 :•  ^ 

1 PhMGE 

I F • M . EQ . 0 ■'  UR  IT  E *:■  - 43n  ' H 1 - P ' 1 ;p . TMP  - DP  » RH  < WH ' L :■  i WO  ■:  L ' » V J S ■ 
IRRMGE 

I r ■:  M . GT . 0 REhD  5 ' 429  ' Z ' L - P ( L > ^ TMP  - DP  ? RH  < UH ' L ' - WO  ■:  L ■' ' 

1 RHRZE ' L 

J=IFIX<Z'.L.^  + 1 . 0E-b)  + t . 

IP  ■ N.  EQ.  0 :'Z<L  ■=-Hl 

I F p;  Z '■  L . GE . 25 . 0 • J= 1 L -25 . 0 ;p  5 . 0 + 2b . 

I F ' Z '■  L ' . GE . 50 . 0 ' J = C Z L.  .j  -50 . 0 ) ' 20 . 0 + 3 1 . 

I F •:  Z < L > . GE . 70 . 0 > J=  f Z L '-70. 0’:>  .-30, 0 + 71. 

IF';.J.GT.33>.J=33 
FhC=Z'::L.:'-FLOhT':  .J-1 
IF',  J.LT.26>  GO  TO  500 
F ri  C = ' Z '■  L . ' - 5 . 0 + F L U H T i - 1 G ) - 2 5 . ) ' 5 . 

I F J . GE . 3 1 F RC=  '■  Z < L ' -50 . O > .-"lO . 

IF  ( J.  GE.  32)  r R C = < Z ( L - 7 0 . 0 > ' 3 0 . 

I F ' , F R L . G T . 1 . 0 F H L =1 . Pi 
.-^  = J+1 

T'  L;>  = TMP  + 273.  15 
TT^273.  15/T<L;.' 

IF'RPI.LE.0.0'  TT  = 273.1 5 / 273 . 1 5 + DP ) 

I F ' WH  < L . LE . 0 . 0 •'  UH  < L =F  ( TT) 

I F ' RH . GT . 0 . 0 ' WH  L ) =0 . 0 1 *PH*WH L > 

IF  ' Z ' L ' . GE. 5 . 0>  GO  TO  520 

I F RHRZE  L ) . EQ . 0 . 0 ' RHZ2  ( L ) =HZ2  '■  J ) * '•  HZ2 P. ) ■'  HZ2  ( J ) > + + FRC 
IF  RHRZE (L  ' . EQ.  0. 0)RHRZE'::L)=HZ1  ■rj>  + <HZl  '■  K ' HZl  ( J> 

IF' MODEL.  EQ.0.:'  GO  TO  S 
IF'  L.  EQ.  1 )WRITE':C,  441  ' 

WR I TE  • 6 - 429  .> Z ( L > > P ( L ' TMP . DP  ^ RH , WH  ■ L • » WO ' L ■ i RHRZE ' L > 

CONTINUE 

IM=0 

NL=HLDRT 

NOTE  THRT  Z-;:L)  MRY  NOT  CORRESPOND  TO  THE  VRLUES  GIVEN  FOR 
MODEL  RTMOSPHERE 
GO  TO  300 

IF ' RRNGE.GT.0.0  ' GO  TO  580 
I F (.  H 2 . G T . 0 . 0 . R N D . H 2 . 1.  T . H 1 ' GO  TO  16 
GO  TO  8 
ITYPE=2 

B E T R = R C 0 S < 0 . 5 + ( F'  H N G E * F'  Fi  N G E ■ ' < X 1 + X 2 ) - X 2 ■ " X 1 - X 1 /'  X 2 > ) C R 

IFCBETR.EQ.O. ) GO  TO  6 

GO  TO  16 

BET=CR#BETR 

X2=RE+H2 

RHGLE=RTRH CXl’-i^S I N '■  BE  T ) (, X2*C0S'l BET ) -IX 1 ) ) ■ CR 
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c 

RhNGE=!s2'!^S I H ' BET  SIN'.  HNGL..E*CR) 

U1S9 

r 

EET^BETfl 

1 80 

c 

GO  TO  3 

0191 

6 

RflNGE*(.X2/Xl  •**2-'.3IN(HNGLE*Cfl.O**2 

0192 

I F k RRNGE . GE . 0 . 0 ' PRNGE  = X 1 * SORT  <. RRHGE  ' -RBS COS < RHGL.E*CH ) ) ) 

0192 

r 

I F ' RNGLE . NE . 0 . . OR . RNGLE  . NE  .180.  ) BET  = RS  I H '■  RRHGE*S I H ■:  RNGLE#CR ) X2 ) 

0194 

IF'  RNGLE.LT.0.  ) HNGLE  = FiNGLE  + P I 

0195 

I F ' RRNGE . LT . 0 . 0 > RRHGE^^-PRNGE 

0196 

bet=bet-cr 

0197 

NR  I TE 6 . 428 ) H 1 - H2  - RflGLE  j RRNGE » BET  < V I S 

0193 

3 

CONTINUE 

0199 

SUriR=0. 

0200 

I F < I XV .LE.2>  PERU ' 5 > 40b > V 1 j V£ , HV 

0201 

I F •:  I KY . LE . 2 ) NR  I TE ' 6 » 406 ) V 1 , V2 » DV 

0202 

I F I TYPE . EG! . 1 NR  1 IE 6 - 407  > HI,  RRNGE 

0202 

I F I TYPE . EQ . ^ NR  I TE  6 , 403 HI,  H2 , RNGLE 

0204 

IF'  I TYPE. EQ. 3)  NRITE'6,409)  HI, RNGLE 

0205 

IF' NOriEL.EQ.0  .'  M=2 

0206 

I F ':  V I S . GT . 0 . 0 :'  NR  I TE  < 6 , 4 1 7 ) V I S 

0207 

I F ( V I S . LT . 2 . 0 . RND . V I S . GT . 0 . 0 ) NR  I TE ( 6 , 442 ) 

0203 

IF'  M.EQ.1'  WRITE’ 6, 414)  N 

0209 

I F ■ I HR2E . EQ . 0 ' NR  1 TE ' 6 , 426 ) 

0210 

I F C V I S . LE . 0 . 0 . RND . I HRZE . GT . 0 ) NR  I TE  ( 6 , 4 1 6 ) I HRZE , HZ ';  I HRZE , L ) - L 

0211 

1=1,2) 

0212 

RVW= 10000. /VI 

0213 

HLRM= 10000. /V2 

0214 

NR  I TE  < 6 , 4 1 8 ) V 1 , V2 , D V , RLRM , RVW 

0215 

RVW=0. 5E-4*':V1+V2> 

0216 

RVW=RVW**2 

0217 

CU=77.46+.459*fiVW 

0216 

CW=43. 487-0. 3473*RVW 

0219 

9 

IF'.JP.EQ.0)  WRITE<6,427) 

0220 

IF(  ITYPE.ECl.  1 ) GO  TO  15 

0221 

no  11  K=1,S 

0222 

VH(K)=0.0 

0223 

1 1 

CONTINUE 

0224 

BETR=0.0 

0225 

SR=0.0 

0226 

IP  = 0 

0227 

c 

0226 

**  NON  DEFINE  CONSTRNT  PRESSURE  PRTH  QURNTITES  EH<l-4) 

0229 

c 

0230 

Y=CR**RNGLE 

0231 

SFHI=SIN(Y,' 

0232 

R1=(RE+H1 )*SPHI 

0233 

IF':Hr.  GT.Z(NL))  GO  TO  13 

0234 

GO  TO  15 

0235 

13 

X=fRE+2<NL>)/  '..RE  + Hl ) 

0236 

IF- SFHI.GT.X:’  GO  TO  14 

0237 

H1=Z<NL) 

0238 

L1=NL 

0239 

3PHI*SPHI'X 

0240 

RNIjLE*  1 80 • 0-RS I N SPH I ) ■ ' CR 

0241 

Pl  = <R£+Hn*SPHI 

0242 

GO  TO  15 

0243 

14 

HMIN*R1-RE 

1 
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0244  UPITE  ' P.  4 3::  ' HMIM 

0245  IPPM=2 

0246  NhME<3'  = 2H5 

0247  CRLL  EXEC' 8<HhME  ■ 

0248  15  CONTINUE 

0249  NOME' 3)  = 2H2 

0250  :hLL  EXECCS-NRMt  ' 

0251  *HLL  MPDfi 

0252  16  UP  I TE'. 6. 445.;' 

0253  STOP 

0254  400  POPriRT  ' 8 13  > F 1 0 . 3 ' 

0255  40 1 FOPMRT ' SE 1 0 . 3 ' 

0256  402  FOPMRT ' F6 . 1 > E 1 0 . : - F6 . 1 ? 2E 1 0 . J > 

0257  404  FOPMRT' 15F5. 2 > 

0258  405  FOPMRT' 8E9. 2' 

0259  406  F OPMRT ' 7F 1 0 . 3 •' 

0260  407  FOPMRT'  '/,28H  HOPlSONfRL  PRTH^  RLT1TUHE  =»F7.3j11H  KM,PRNGE  =- 

0261  1F7.3«;3H  KM) 

0262  408  FOPMRT'. 'SCRNT  PRTH  BETWEEN  RLT1TUDE3  HI  RND  H2  WHERE"- 

0263  l."Hl  ="..F7.3^8H  Ml  H7  =.F7.3MSH  KM, ZENITH  RNGLE  =,F7.3, 

0264  2SH  DEGREES  :' 

0265  409  FOPMRT':: M 39HSLRNT  PhTH  TO  SPRCE  FROM  RLTITUHE  HI  =,F7.3- 

0266  r'KM",/,  "ZENITH  RNGLE  ^",F7.3,SH  DEGREES;:' 

0267  414  FOPMRT': /1 0X- 18H  MODEL  RTMOSPHEPE  MM21H  = 1962  US  STRNDRPD  ,' 

0268  416  FORMAT':/ 10K,  ISH  HRZE  MODEL  , I 1 , 3H  = ,2R2,13H  VISUAL  RANGE  :' 

0269  417  FORMAT':. ./25X"HR2E  MODEL  =",F5.1,"  KM  VISUAL  RANGE  AT  SER  LEVEL"',) 

0270  418  FORMATS/ 1.20HFREQUENCY  RANGE  Vl=  ,F7.1,13H  CM-l  TO  V2=  ,F7.1, 

0271  15H  CM-l ,./,  "FOR  DV  =",F6.1,8H  CM-l  '::,F6,2,"  - ",  F5. 2,  "MICRONS) " ) 

0272  426  FORMAT':.- laX,  "AEROSOL  SCATTERING  NOT  COMPUTED,  IHAZE-0"  ' 

0273  427  F0RMAT':/.'//./9./31X,  20H  HORIZONTAL  PROFILES.-') 

0274  428  FORMAT '■/.-•  "HI  = ",F7. 3,  "KM,  H2=" , F7. 3,  "KM,  ANGLE= " , F8 . 4 , " , GEOM.  R" 

0275  r ANGE  =",F7.2,  "KM".-  "BETA=" , F8. 5,  "DEG.  , VIS  = ",F6.  1 - ) 

0276  429  FORMAT':3F10.3,2F5.  1,2E  10.3, 2F1 0.3) 

0277  430  FORMAT'  10X,  " INPUT  METEOROLOGICAL  DATAf".- , " Z = ",F7.2,"  KM,  F-",F7 

0278  1.2,"  MB,T",F5.1,"  C,  DEW  PT . TEMP " , F5 . 1 , " C,  PEL  HUMIDITY=" , F5. 1 , 

0279  2 " •. ' " H20  DENS  I T'i'= " , 1 PE9 . 2 , " GM  M-3 , OZONE  DENS  I TY  = " , E9 . 2 , " G 

0280  3M-  3 " - " V I SURL  RANGE= " , 0PF6 . 1 , " KM , RANGE= " , F 1 0 . 3 " KM  " !' 

0281  431  FORMAT' 4':;F6.2,2F7.5)  ' 

0282  433  FORMAT(;/.--"TRAJECTORY  MISSES  EARTHS  ATMOSPHERE.  "CLOSEST  DIS" 

0283  r’TAHCE  OF  APPROACH  IS" , F10. 2,  IX,  - , IX,  "END  OF  CALCULATION":' 

0284  441  FORMAT'- /./"MODEL  AIMOSPHERE  NO.  2",/  4X,"Z  'IKM) " , 3X,  "P  ':MB;)",4X, 

0285  1 T (C:>  DEW  PT  XRH  H2Cii:GM.M-3)  03'::GM.  M-3)  NO.  DEN.') 

0286  442  FORMAT'  -/"FOG  CONDITIONS  MAY  EXIST  AT  SEA  LEVEL  FOR  THIS  VISUAL  Ri 

0287  IN'UE",/,"  IF  SO  THEN  ASSUME  THE  TRANSMITTANCE  DUE  TO  FOG  IS  GIVEN 

0288  SBr"'.-"  THE  TRANSMITTANCE  AT  0.55  MICRONS":.' 

0239  425  FORMAT*.  10X,7F10.3) 

0290  445  FORMAT'  //"#**ILLEGRL  INPUT  PARAMETERS  FOR  NRDA*#*"..- 

0291  1"+*4CALCULATI0NS  TERMINATED***") 

0292  END 


**  NO  ERRORS**  PROGRAM  = 0285 1 COMMON  = 05533 


A -7 


I 


f-  Mljt 


KTN4  lump  1 1 EP:  HP24177 


^ 'SEPT.  1974) 


UUU  1 

nn.L 

uuu2 

PPULPRM  self:  '5'90;' 

Ut10  I' 

III  MENS  I ON  NRMF.'LS  ' 

0004 

iJiMMON  TEMP'-  40  ' . PPES'  4ii  ' - '.'1  - V2l  DV-  IPhDj  1 Mn:7-  RHZ2':  20  ' 

000'=^ 

i:  iiMMON  PRNGE  - NU'  4n»  8 

000  L 

L OMMON  Z 34  ' ' P '■  34  ' . T '.  .;4  ' < E H '.  8 < 34  :'  < NH  ('  34  .'  - M ? Nl  j PE  ' i U • L 

0007 

1 i iMMON  no  ^54  ' - HZ  1 '.  34  ' < HZc  ''  6 ^ TX  (10:'^  VH  ■ S ' ’ W 8 > - E 7 8 ' 

000  c 

COMMON  C4 ' 1 33  ' > 05  • 1 5 ' ' RHRZE  ■;  34)  ? VX  • 45  • - C7>:  45  ' 

0009 

C OMMON  C7R ' 45  . , L I'P  '■  40 - Rl..  T '■  40 .'  ^ I PPM 

0010 

L UMMON  X 1 - X2 . H 1 < H2 . N ^ NF'  ^ HP  1 - NP2  ? TX  1 ^ T 7:2 , )'N  < VH  1 , ) N2 

001  1 

L UMMON  1 . Ip . L 1 , L2  ^ PS  ^ T S < X - PEP  - I P ^BPl  E V - T MP  - P2 . H ' US 

00  1 2 

L OMMON  PS  I - SF'H  1 , 1 RET  - THETR  < PH  I < BE  T . BEl  H , SRLP  i SP 

00  1 3 

LijMMuN  HM  I N ) L 1 ’ HM ' P:V  ' L Pi ' PN  n LEN » SUMR  ? RB ' RLRM 

0014 

L OMMON  I X'l' . V I S - I V 1 - I V2  - I HRZE  - I TYPE  < RNGLE  ’ F I LL ':  1391  ' 

0015 

r OMMON  MM , MOHEL  < T 1 - T2  ^ TM I N i P 1 - P2 ' PM  I N 

0016 

DRTR  NRME- 2H0E<  2HGT, 3H3 

001  7 

GO  TO.  15- 19'  I PPM 

00  1 S 

15 

LiO  17  L=1-NL 

00  1 9 

P-=P'  L '.'1013. 0 

0020 

TS  = 273.  15  T'-:i  ' 

002  1 

X=^PS>TS 

0022 

pt=ps*supt<;ts) 

0023 

li  = 0.  1+UH'  L ' 

0024 

EH'  i-l::'=0.s*pt>:: 

0025 

PPH  = 4. 56E-5*Ii*T':  L ' 

002  b 

EH ' 2 * L ' = Ii4pPU'i*-E7:P ' 6 . 08* ''  29t-  .0  T > L ' ~ 1 . 0 ' ' 

0027 

EH  ■ 8 ’ L ) = D*  <:  PPU  + C1 . i 2^^  ' PS-PPW  ' ■ ♦ E XF’ ' 4 . 56  ♦ ' . ti  T '.  L ' - 1 . 

0023 

EH'  5»L  ' = . 125E-2*.:-HH'  1 ' 

0029 

EH':6'L ' = .467E-3*::*H0  1 ' 

0030 

EH'-  3»L  ;.=X 

0031 

HRZE=HZ1 CL) 

0032 

IF'M.EQ.2)  HRZ£=hHhZE  1 ■ 

0033 

IF' Z':  1.  ' . GE.  5. 0 J lO  TU  1=10 

00  34 

IF' M.EQ.2.RND.  IHhZE.EO..  • HM..E=h.  . i.  . 

0035 

IF'  IHRZE.EQ.Z.RNIi.M.Er..  HHZE  = iHMZ: 

0036 

I F ' V I S . LE . 0 . 0 ' GO  T 0 1 m;!  , 

0037 
00  3S 

IF' M.  EO.  2)HRZE=  r . SB'-*'  ' ''  '"'I  - '*^R  1 '-L  ' 5.0 

IF  ' M.NE.2.<  GO  Tu  15u  ' 

0039 

0040 

150 

HRZE=6. 389*'  HHZ2' 1 ' uhh  -^^,  "i  ^hHm.'E  i ' 5.0-HHZ2C 

IF'.  HRZE.  LT.  0.0  ' HhZE  n.u 

0041 

EH'  4>  L)  = 3. 53c6E“4*HR.’E 

0042 

IF' M. EQ. 2 ' EH' 4. L '=HRZP  hhRZE • 1 ■ 

0043 

E H I-  7 » L ) = 1 . ti 

0044 

FEF-1 . 0E-6^CCO^X*  101  ,:.0  273 . 1 5-F  PH-CN 

0045 

IF  a.EQ.NL)  GO  TO  16 

0046 

IP' MODEL. EQ. 0. RND. L . GE. 1 ' GO  Tu  26 

004T’ 

T2  = TCL+1  ' 

0048 

H2=WH<L+1  .' 

0049 

PPH=4.56E-6*U2^T2 

0050 

EH '.  7-  L ' =0 . 5* ' REF  + 1 . 0E-6*- ' CO*F' 'I L.+ 1 :■  .■-T2-PPn*i:  W ) ) 

0051 

16 

IF  CL.EO.HL)  EH'I7-L)=0. 

0052 

IF  'Hl.GE.ZCL,:')  L1=L 

0053 

I F ' .JP . EQ . 0 ) I'JP  I TE  ( 6 !.  434  L > Z C L i"  EH C K - L ) - K = 1 > 8 ) » REF 

00154 

EH.:  7,L)=EH(7-L  ) + l .0 

0055 

17 

CONTINUE 

0056 

170 

IP  = -1 

I.  n p 1 


li  ■'  > 


-hZ2'.L  ■ 2* . U 
L 23 . 0 ' 


A -8 


J 


I-  RGE  yutiE’ 


RR  57 
Ry58 
y058 
8080 
tujb  1 
tiut'2 
0087-! 

0084 

0085 
0068 

0067 

0068 

0069 

0070 

0071 

0072 
007  3 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 
0033 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 
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II  =0 
X1=H1 

ChLL  point  ' hi  - VN«  L,  NPl  , T7MP> 

11  = 

9 1 t:7'10. 

L1=L 

t;:i  = ix<7.-' 

DO  18  F=l«8 

18  E' f ;'  = t:-:'.p  ■ 

LBR  = 0 

IF  ITVPE.  EO.  1 ■ GO  TO  Gb 
IF  ITYPE.EQ.  3:.'  H2=Z'  HL) 

I F ' HNGLE  . GT  . 90 . 0 GO  T 0 28 

19  IF  ' RNGLE.  GT.  90. 0.  RHIi.  NPl . GT,  0.'  Ll=LUi 

12  = NL 

IF  ' ITi'PE.Eu.3;'  GO  TO  20 
ORLL  POINT  ■:  H2 . YN . L n NP - TX ^ I P 
T2  = TX^.9> 

92  = TX':'10;' 

L2=L 

IF  '-.NP.GT.O;'  L2  = L2-1 
c H ' 8 ^ L 1 .j  = E ' . 8 !' 

20  no  21  K=1.6 

EH  ■:  K - L 1 :■  =E  • K ' 

IF  'IT Y P E . E Q . 3 .j  1 1 0 T iJ  2 1 
EH'IK.L2-t-l>  = TX^K::' 

21  continue 

I F C I TYPE . HE  . 3 > EH  ■ 8 . L2  1 .)  = TX  ( 8> 

I'-  ai.EQ.L2'  TX1=TXUYN~EHC7»L1  :' 

0 

,:♦*,**  MOW  DEFINE  VERTICRL  PRTH  QUANTITIES  VH'l-6.:' 

C 

IF  'i  JP.EQ.O)  WRITE'IEj  420  j 

DO  25  L=L1«L2 

X1=Z<L> 

X2=Za-H  > 

IF  '::L.EQ.L1'  X1=H1 
IF  a.EQ.LS'  X2=H2 
DZ=X2-X1 

IF  (L.EQ.NL)  DZ=Z(L>-Z':;L-1  > 

Ds=riz 

c 

UPWARD  TRAJECTORY 
C 

Ri RE  + X 1 ''  ( RE-*-Xl2 
THETA=ASINC3PHI  I'.-CA 
PHI  = ASIN'::SPHI*RX..' ■CA 
BET=THETA-PHI 
SALP=RX*SPHI 

IF  (SPHI.GT.  l.E-10>  DS='IRE-hX2>*SIH(BET*CA>/SPHI 

BETA=BETA-i-BET 

PSI=BETA-t-PHI-ANGLE 

PHI=180. -PHI 

SR=SR-t-DS 

LL  = L-LULBR-H 

DO  244  K=l»8 
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y 1 1 

0114 

0115 
0 1 1 
0117 
01  li; 
0 1 1 4 

0 1 2y 
0121 
0122 

0123 

0124 

0125 
01 2'? 
0127 
012i: 
0129 
01  20 

0131 

0132 
0132 

0134 

0135 
013c. 
0137 
0 1 3S 

0139 

0140 

0141 

0142 

0143 

0144 

0145 
014b 
0147 
0143 

0149 

0150 

0151 

0152 

0153 

0154 
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EV  = DS4EH':.K,  L ' 

IF  'L.EQ.NL'  GO  TO  22 

IF  EH ' K ' L ■ EO . 0 . 0 . OP . EH K ? L+ 1 . EQ . 0 . 0 GO  T 0 23 

IF  's E H ' , K ’ L > . E I! . EH''  1'  » 1.  + 1 ' ,'  0 U I 0 '2 4 

E = D b E H ' 1^  ' 1 . ' - E H ' f > L.  h 1 ■'  ' ' H L Li  G ( E H < f . < L. E H K j L>  1 ■' 

GO  TO  24 

22  IF  ■ EH • K • L . E 0 . 0 . 0 ' GO  7 0 23 
IF  EH' K' L-1  ' . EO.  0. 0 ' GO  TO  23 

IF  . EH '.  F ' L ' . EL' . EH ' f'  ? L - 1 ) ' OO  T 0 24 
EV  = EV  HLOG' EH' F,L-1  ■'  EH'::F-I.  > ' 

GO  TO  24 
2:;  EV=0. 

24  VH' K >=VH' K >4EV 
f 

MpWRPIi  PHTH---STOPE  PhPHMETEPS 

r 

244  UN'  LL-f  .'  = EV 
^VP'iLL'  = L 
HLTCLL'  = XI 

TEMP  ai  ' ^ SORT ':  T L ) *T  ( L+  1 ) 

ppEsai ' = soPTc p l > *p '.  l+  i ;■  > 

I F ■:  .IP  . EQ . 0 ;'  HP  I TE '"  6 ' 435  7 L . X 1 < i„  VH  ■ K 7 - K=  1 . b - BETH  ^ THE  7 H . '3P 

IF  'L.GE.NL  ' GO  10  25 

IF  ' L+  1 . EQ . L2  ' EH ' 7 j I.  + I ) =VN 

IF  vL.EQ.Ll  :.  EH'.  7iL)=TXl 

PN  = EH':.7<L+1  ' EH'X^a':' 

SPHI=SPHI*RX  PM 
IF  'ISHLP.  GE.  PM;'  ■BPHI'^SRLP 

25  COOT  I HUE 

I PPM  = 1 

GO  TO  100 

2b  IPPM  = 2 

GO  TO  100 

23  I PPM  = 3 

100  NRME'3'  = 2H.: 

CRLL  E'-'EL  (3- HRME  ' 

CRLL  MPDR 

420  FOPriRT  . • ■••'25X»  "VERT  ICflL  PROFILES" » -'59X»  "BETH  THETR" 

13!!'  " PRHGE  " ' 

434  FOPflRT'-;  I3jF5.  1 . 9'! E8. 2)  :> 

435  FOPMRT  ' 1 3 ' F5 . 1 , EES . 2 ? 2FS . 3 j F6 . 1 7 
EMU 


COM NON 


055 


HiJ  EPP0Pb4  4 


PPOGRRM  - 01831 


0001 


MH4  compiler:  HPE4177  '..GEPT.  l'?74> 


7i00l  FTfM 

0 0 0 G P R 0 G R H t1  G E i j T 7 5 « 0 ' 

Oooj  rilMENSIOH  NRMl:  (' 3 ' 

0004  r uMMOM  TEMP  ' 40  ' < PRES ' 40 « V 1 1 V2 ' D V . i RhDH  < L I'IhX  « RHZ2  ■:  20  ■ 

0005  rriMMON  PRUGE- NW' 40' 8 ■ 

oGOt'  L-ijMHijH  Z ' o4  ' P c .34  ' ' T ' '4  < EH  8 34  :i  IJH  >'  J4  ^ H ' ML  ' PE  ? UN  ’ L'L  < P I 

000  1 L I iRI'lLiH  NU  V 34  ' ’ HZ  1 ' >4  ' ’ HZ3  8 ) ' TX 1 O'  ' » VH ' 8 ' < W 8 ' ’ E 8 

0 0 0 P C I j M M 0 R C 4 1 .3  3 < C 5 ' 15  ’ ' hi  H R Z E ( 3 4 » V ':■< ' 4 5 ' C 7 ' 4 5 > 

0009  COMMON  C7R  ■.  45 < L VP  ■:  4n  > - RL  T ' 40  > ? I PPM 

0010  LUMMON  XI  , X2  - HI  - H2 ' N - NF' - NP 1 7 NP2-  TXl  7 TFIC  7 YN  7 VN 1 - VN2 

00 1 1 COMMON  L 7 JP  7 L 1 7 L 2 ^ P3 < TG  7 X 7 REF  7 I P 7 LBR  7 EV  7 IMP  7 K2  7 H , DS 

00 1 2 COMMON  PS  I 7 9PH I 7 TRET  7 THETR  7 PH  I 7 BET  7 BETR  7 SRLP  7 SR 

00 1 3 C OMMON  HM I N 7 l.L  7 HM 7 RX  7 CR  7 PN  7 LEN  7 SUMR  7 RB  7 RLRM 

0014  I.  OMMON  I XY  7 V'  I S 7 I V 1 7 I Yj , I HRZE  7 I TYPE  7 RNGL  E 7 F I LL  <1  39 1 :■ 

00 1 ^ COMMON  MM  7 MODEL  7 M 7 Tx • TM I N 7 P 1 7 P2  7 PM  I N 

0018  DRTR  NRME  ■ 2H8E 7 ERG T 7 8H4 

0 0 1 I ' j 0 T 0 4'  5 7 487  4 8 7 3'  5 * I F'  F'  M 

0018  25  LBR  = L2-L1+LBR+1 

0019  I PPM  = 4 

0020  L 9L-  EXEC  3 7 NRME  ' 

0021  C 

0022  ROPIZONTRL  FRTH 

0023  C 

0024  28  Di.i  27  K=l7  8 

0025  W 1 ■=RRNGE*EH'.  K<  1 ' 

0028  IF  - MODEL.  GT.  0"'  N' K..'  = RRNGE+TX(F:;) 

0027  27  UN^  1 .K  ' = Wa-  ' 

0028  LCR'l'  = LI 

0029  IF' MODEL.  EQ.0;'  TEMPO:'  = T(l> 

00  30  IF' MODEL.  NE.  O ' TEMPO  ' = T1 

00  31  IF  MODEL . EQ . O ' PRES 1'  = P ':  1 .' 

0032  I F ' MODEL  . HE  . 0 PRES  O'  = PI 

0033  1 1 ’■  MODEL  . EQ . 0 RL  T O ' Z 1 ) 

0034  IF' MODEL.  HE.  0:'  RLT'l'  = HI 

0035  LBP  = 1 

0036  GO  TO  49 

0037  28  CONTINUE 

0038  C 

0039  DOHNURRD  TRRJEC 1 ORG' 

0040  C 

0041  h.2  = 0 

0042  IF  'NPl.EQ.l'  L1=L1-1 

0043  L2=L1+1 

0044  YN1=YN 

0045  L0=L1+1 

0048  IF  ' H2.  GT  . Z':  L 1 + 1 , OR.  HI . EQ.  H2>  GO  TO  30 

004  7 f F ■ HP  1 . EQ . 1 . RND . H2 . GE . Z L 1 + 1)  > GO  TO  30 

0048  CRL  L PO I HT  ■!  H2  7 YN  7 L 7 NP2  7 TX  7 I P > 

004  9 T2  = TX':.9> 

0050  P2  =■  rXO0> 

0051  DO  29  K=l78 

0052  29  H' F )=TX':.K.' 

0053  TX2=TX'i7:> 

0054  i'H2=YH 

0055  IF  <H2.LT.Hl:'  H=H2 

0056  L2=L 
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0057 

IP  ai.EQ.L2'  Tj:2^TXl+'.'N2~EtK7,L) 

005? 

IF  '.H2.GT.H1'  TX1  = IX2 

0059 

IF  ' L1.EQ.L2.HND.H2.L1 .HI ' YN1-TX2 

0060 

30 

R0=':  PE^Hl  >*SPHI*VN1 

0061 

IF  ‘ H 2 . G E . H 1 ' Y N 2 = 1' N 1 

0062 

no  31  L=1.L1 

0063 

HrilN=H0 'EH' 7- 1 >-Ff 

0064 

IF  ' L.EQ.Ll  ' HniN^H0  YM1  -RE 

0065 

LM1N=L 

0066 

IF  ■ HN!N.  L£.  l^  L-^  1 ‘ • Gu  TO  32 

00157 

31 

CONTINUE 

0068 

id. 

X = HMIN 

0069 

IF  'HNIN.LE.O'  GO  TO  34 

0070 

LHLL  POINT  (X? YN1L<NP<TX<  IP) 

0071 

TMIN  = TK':9) 

0072 

PM  IN  = TKM0  .' 

0073 

LMIN  = L 

0074 

TX3=T>:‘  7 :■ 

0075 

IF  L2.  EQ.  L.  OP.  [.  1 . EQ.  L ' T X3=YN2+TX < 7 :■ -EH 7 < L 

0076 

IFaX3.LT.0.O:'TX3  = TX'  7 ' 

0077 

IF  ' L1.EQ.L.RND.H2.GE.H1  ' GO  TO  33 

0078 

HMIN  = R0  - TX3-PE 

0079 

I-  ' RES' X-HMIN > . GT . 0. UU01  ' GO  TO  32 

0080 

35 

IF  •; L 1 . EQ . L . RND . H2  . GF  .HI'  YN 1 =TX3 

0081 

IF  (L2.EQ.L.RND.Ll,NE.l.2>  YN2  = TX3 

0062 

IF  ' H2.GE.H1  ' TX2=TX3 

0083 

IF  ':H2.GE.H1  ' L2=L 

0084 

IF  <H2.GE.Hl.OR.H2.LI.HMIN.'  H=HMIN 

0085 

l'JRITEf6<  436  .'  HMIN 

0086 

IF'  H2.LT.HMIN,.'L2=L 

0087 

IF  ' H2.  LT.  HMIN  ' HP  I TE ':  6 < 440  ) HMIN 

0088 

GO  TO  35 

0089 

34 

MRITE' 6< 436)  HMIN 

0090 

IF  '■  H2.  LT.  HI GO  TO  35 

0091 

IF  aTYPE.EQ.3.0P.H2.GE.Hl)  HP  I TE ':  6<  437 ) 

0092 

IT'i'PE  = 2 

0093 

TY2=EH(7<  1 ,' 

0094 

LMIN=0 

0095 

L2=l 

0096 

H2=0. 0 

0o97 

H = 0. 0 

0098 

c 

0099 

NOW  DEFINE  VEPTICRL  PRTH  QURNTITIES  VH(l-4) 

0100 

c 

0101 

35 

IF  UP.EQ.0)  UPITE<6<420) 

0102 

L = L0 

0103 

LL  = LBP 

0104 

DO  40  1=1 <NL 

0105 

LL  = LL+1 

0106 

L=L-1 

0107 

REF=EH'  7<L;' 

0108 

IF  I.EQ.l)  Pr.r=-iNl 

0109 

IF  <I.EQ. 1.RND.K2.EQ.  1 ' REF=YN2 

0110 

IF  ' L. EQ. L2. RND. K2. EQ. 0 ’ REF=TX2 

011  1 

IF  ':l.NE.l)  Xl=2a  + 1' 

0112 

X2=Z(L:' 

A-12 


0 


y 1 1 1; 

yl  14 
yl  15 
y 1 1 fc 
yi  17 
y 1 1 £. 
y 1 19 
yi20 
y 1 2 1 
0122 

0123 

0124 

0125 
y 1 2 b 

0127 

0128 

0129 

0130 

0131 

0132 

0132 

0134 

0135 
013t‘ 
0 1 37 
0138 

0133 

0140 

0141 

0142 

0143 

0144 

0145 
0148 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 
015b 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 
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IF  a.EG.L2.HND.F2.EQ.0>  X2=H 

IF  ' L.EQ.LMIM.HNLI.K2.EU.  1 ' X2  = Hf1IN 

HM=' PE+Xl >*SPHI-RE 

IF  ':  HM . GT . 2 ':  L ' . AND . HM . GT . X2 X2=HM 

RX  = '-  RE  + Xl  ' -'  RE+XS  ' 

D3=Xl-::2 
HL-P  = 90. 0 

THET=hSIN'..SPHI  ' ■ Cm 

shlp=p::*sphi 

IF  ':  ABSi^X2-HM  ‘ . GT,  l.yE--5  ' hLP  = hSIN'  ShLP  ' CR 
BET=hLP-THET 

IF  ( SPHI . GT.  1 . yE-10:'  D8=^  RE+X2;'+SIN''BET*CR)/SPHI 

THETR=1S0.0-THET 

BETR=BETH+BET 

P3I=BErR-RLP-HNGL  E + l;;;0. 0 

3F-SR+DS 

DO  399  K=l?8 

R L = E H '%  K * L ■' 

BL  = EHO:iL  + l ' 

IF  'L.EQ.Ll  ' BL=-E'K> 

I>^  ' L.  EQ.  L2.  RHD.  H2  . L I . HI . RHD.  H2,  GT  . U.  O ' RL  = W';  K';' 

IF  a.EQ.LMIN.RND.H2.GE.Hl>  RL  = TX':K;' 

IF  a.EQ.LMIH.RND.REi:'.H2-HM;'.LT.  1.0E-5'  RL  = TX(K> 

'K2.EQ.0'  GO  TO  36 
IF  a.EQ.L2;'  BL  = N'K> 

IF  a.EQ.LMIH,:'  RL  = TX' K • 

36  IF  f RL . EG . 0 . 0 . OR  . BL . Eu . 0 . ' GO  TO  38 

IF  (RL.EQ.BL::'  GO  TO  3 7 

EV  = ri3* '.  RL-BL  ■•' /RLOG  '■  RL  -BL  > 

GiJ  TO  39 

37  EV=D3*RL 
G*'i  TO  39 

38  E'.‘--0.0 

39  VH'-  K ;'  = VH'  K;'  + EV 

C 

DOHNWRRD  PATH  — -STORE  PARAMETERS 
C 

39-  WIML,F'  = EV 

LBR  = LL 
L'-’R'LLi  = L 
RlT':LL'  = XI 

TEMPED  = SORTvTa>  + n:L+l  >> 

PRES  aL  > = SORT ' P ' L ) ^P ' L+ 1 > > 

IF  ^JP.EQ.0)  NR  HE  <6 '435;'  Lj  X1 1 (VH'::K '»  K=1 , 6>  > BETA,  THETA,  SR 
IF  a.EQ.L2.RND.H2.GE.Hl>  GO  TO  45 
IF  a.EQ.LMIN.RND,K2.EQ.  1.)  GO  TO  43 
IF  ':.L.NE.1>  RN»=REF  EH' 7,L'l) 

IF  <L.EQ.L2+1>  RN=REF.  TX2 

IF  a.EQ.L2.RHD.F2.EQ.0>  RN=REF/YN2 

IF  a.EQ.  (LNIH+l  '.RND.K2.EQ.  1>  RN=REF- TX3 

IF  '^SRLP.GE.RH)  RM=l.y 

SPHI=SRLP*PN 

IF  ':.L.  EQ.  L2.  RHD.  K2.  EO.  0>  GO  TO  41 

40  com  I HUE 

41  I PPM  = 1 

TEMP'LL)  = SQRTa2*T':L)> 
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0169 

F RES 'LL'  = Si‘'PT'.  R2-P' L ' ) 

0170 

GO  TO  100 

0171 

43 

IPRM  = 2 

0172 

TEMP'.  LL)  = SORT'  f MIH*T'  L+n) 

0173 

RPES'.LL  .'  = SORT' RMIN^P' L+ 1 > ) 

0174 

GO  TO  100 

017? 

4? 

IPRN  = 3 

0176 

TEMP'LL'  = SoRT'LTUT' L' 

0177 

RRES'-LL  ' = S0RT':  P2-^P'.  L ') 

017S 

GO  TO  100 

0179 

49 

IPRM  = 5 

0180 

1 00 

1 hLL  E!sEL-<y!i  NOME  ' 

0181 

ChLL  MPDR 

01 1 82 

420 

FORMAT  /.•'2?JN  "VERTICAL  PROFILES"  »/ 

?9X » 

"BETA  THETA" 

0183 

1 3! ! > "RANGE " .' 

0184 

4 3? 

F 0 R M H T 'L I S' ' F ?' . 1 ; 6 E o ■ S ’ d F 0 . S i F"  6 . 1 > 

018? 

436 

FORMAT  (:.-.'"HMIN  = "jFlO.S) 

0186 

437 

FORMAT  -'"PATH  INTERSECTS  EARTH  - PATH 

CHANGED  TO  TYPE  2 WITH  H 

0187 

r 2 = 0.0  KM"  ' 

0188 

440 

FORMAT  (.■/"H2  WAS  SET  LESS  THAN  HMIN  AND 

HAS 

BEEN  RESET"  ■•■"EQUAL 

0189 

I'TO  HMIN  I.E.  H2  = "-F10.3) 

0190 

END 

i 

iJO  ERRORS--  PROGRRM  = 01863  CONMOH  = 0?5 

PHLiE  C.10O1 
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0U01 

FTIM 

0002 

PPDGP'HM  SEGT 4 ' 5 ' 9y  ' 

0003 

D I MENS  I OH  HRMEc.3  ' 

0004 

COMMON  TEMP ( 40;- , PPES ' 40  i V 1 , V2 ? DV > I RhD i.  LMhX » HHZ2 ( 20 ) 

0005 

COMMOH  RflNGE-l'JWUOiS  ' 

000P: 

COMMON  2('34)»P ( 34 > - T '■  34 > « EH  1 8 » 34 ) . WH < 34 ) - M » HL , PE » CN , C 

0007 

COMMON  UO  ( 34  ' ' HZ  1 ( 34  ' . HZ2  ( 6 ) > TX  a 0 ' - VH 8 ' > W 3 ) < E ( 8 ) 

000S 

COMMON  C4c:133>,C5'.  15  ' - RHRZE<34> , VX-45  ^ C7(45> 

0003 

C OMMON  C7R(45  ' p LVR':.  40)  - RL  T>:  40)  p IPRM 

0010 

COMMON  X 1 p X2  p H 1 p H2  p N ' NP  p NP 1 , NP2  p TX 1 p TX2  p VN - YN 1 p YH2 

001  1 

COMMON  L p JP  p L 1 P L2  P PS ' TS  p X - PEP  p I P p LBR  p E V p TMP  p K2  p H p DS 

0012 

COMMON  PSI p SPHI p THET' THETh- PHI P BETp  BETRp  SRLPp  SR 

0013 

COMMON  HM I N p LL  p HM  p RX  p CR  p RN  p LEH  p SUMR  p RB  p RLRM 

0014 

COMMON  I XY  p V 1 8 P I V 1 p I V2  p I HRZE  p I TYPE  p ANGLE  p F I LL  (.  1 392  > 

0015 

COMMON  T 1 p T2  p TM I H P F 1 p P2  p PM  I H 

0016 

DRTR  NAME  ••  2HSE  p 2HGT  p 2H5  •• 

0017 

GO  TO  '41ji43p45p47p49>  IPRM 

00  IS 

41 

IF  (HMIN.LE.O)  GO  TO  47 

0019 

IF  ( L E N . E Q . 0 ' W P 1 1 E C 6 p 4 3 3 ) 

0020 

IF  CLEN.EQ.O'  GO  TO  47 

0021 

IF  (LEN.EQ. 1 ' NPITE<6p439) 

0022 

I<.2=1 

0023 

Y1=X2 

0024 

IF  <RBS<X1-HMIN '.LE.  0.001  GO  TO  47 

0025 

H=HMIN 

0026 

L-l.2+1 

0027 

IF  (NP2.EQ.r-  L = L-1 

0023 

B=BETfi 

0029 

PH=1 80. 0-RSIN(SPHI >/rR 

0030 

TS=SR 

0031 

P3=PSI 

0032 

DO  42  K=l»8 

0033 

42 

E(K)=VH(K) 

0034 

LSTORE  = LBR 

0035 

GO  TO  35 

0036 

43 

BETR=2.*BETR-B 

0037 

PSI=2.^PSI-PS 

0038 

SR=2.^SR-TS 

0039 

c 

l.ONG  PATH  TAKEN 

0040 

PH1=PH 

0041 

DO  44  L=li8 

0042 

44 

VH  - K.-=2.  *VH(K)-E-  K) 

0043 

c 

0044 

DOWNWftRD  H2.GT.H1--L0HG  PRTH  STORPiGE 

0045 

c 

0046 

LLMIN  = LBR+1 

0047 

LBR  = 2*LBR-LST0PE 

0043 

DO  446  LL^LLMINpLBR 

0049 

LMRP  = LBR  - LL  + LSTORE 

0000 

RLT(LL)  = flL^<LMHP+2:- 

0051 

IFCLL.EQ. LLMIN)  Gu  TO  442 

0052 

TEMPaL.  ^ SQRT(T(LMRP+1  )*T<LMRP  + 2>) 

0053 

PPES<LL)  = SQRKP(LMRP+1  >*P<LMRP'r2)> 

0054 

GO  TO  444 

0055 

442 

RLT(LL)  * HMIN 

0056 

PRES(LL)  = SQPT<PMIN*PaMRP+2>> 

A-15 


:i5pi 
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0057 

TEMP'LL'.'  = SORT' TMIN^T' LMmP  + 2;o 

0058 

444 

CONTINUE 

0059 

DO  446  K = 1 ? 6 

0060 

44b 

NW'  LL.K)  = Wl'KLMRP+1 -P  ' 

0061 

GO  TO  47 

0062 

45 

DO  46  K=1'8 

0068 

46 

VH';  K>=2.  O^tVH'  I ' 

0064 

c 

0065 

DOWNWRRD  H1.LT.H2 Hl.NE.HMIN 

0066 

c 

0067 

LLMIN  = LBR+1 

0068 

LBP  = 2*LBR 

0069 

Du  468  LL=LLNIN,LBR 

0070 

LMRP  = LBR-LL 

0071 

RLTaL.^  = RLT';:lMRP+2' 

0072 

IFaL.EQ.  LLMIN)  GO  TO  464 

0073 

TEMPCLL)  = SQRTa' LMRP+1  'i:T'::LMRP  + 

0074 

PRES LL  ' = Sort  ' P '■  LMRP+ 1 ) *P  ( LMRP  + 

0075 

GO  TO  466 

0076 

464 

RLTaL'  = HMIN 

0077 

TEMP  ■:  LL  ) = SORT  TM I N-T  <.  LMRP+2 ) 

0078 

PRES(LL)  = SORT ( PM  I N^P' LMRP+2)) 

0079 

466 

LVR<LL'  = LYR<LMRP+1) 

0060 

DO  468  K=1p6 

0081 

468 

WW(LL>K)  = WUaMRP+ljK) 

0082 

BETR=2.0*BETR 

0008 

SR=2. 0+SR 

0084 

IF  (H2.EQ.Hi:'  GO  TO  47 

0085 

RN=TX1  YNl 

0086 

SPHI  = SIN(RNGL  E+CR  ' 

0087 

IF  (3PHI  .LT.PN;-  3PHI=SPHI.-  RH 

0088 

GO  TO  19 

0089 

47 

CONTINUE 

0090 

WRITE<6> 406'  HM 

0091 

DO  48  K=1p8 

0092 

N'  P ;=VHO() 

0098 

48 

CONTINUE 

0094 

49 

WRITE  ''6-419' 

0095 

WR I TE  ( 6 > 42 1 ) ' W O - k=  1 j 4 ' - W < 8 ) 

0096 

1 = 1 

0097 

L=1 

0098 

IVl=Vl/5.0 

0099 

IF':  V2.GT.  32767.  ) V2=32767. 

0100 

IV2=V2/5.+.99S 

0101 

rvi=5*ivi 

0102 

I'v'2=5*IV2 

0183 

IF  (IV1.LT.850)  IV1=350 

0104 

IF  (DV.LT.5.  ' DV=5. 

0105 

HRMF-fS)  = 2H5 

0106 

GO  TO  100 

0107 

35 

IPRM  = 4 

0108 

NRME(8)  = 2H8 

0109 

GO  TO  100 

01  10 

19 

IPRM  = 2 

0111 

NHME'iS)  = 2H2 

0112 

100 

CRLL  EXEC(8>NRME' 

A-16 
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0113 

0114  40b 

0115  419 
0 1 1 

0117 

0113  421 

0119  433 

0120 

0121  439 

0122 

0123  500 

0124 


ChLL  MPDR 
FOPMRT'.7F10. 3) 

FORMRT  (.lbX»3SH  ECU  I LRVENT  3ER  LEVEL  RBSORBEP  RMOUNTS/- 
17X"  NITROGEN  'rCONT'  H20  (CONT)  MOL  SCRT  REROSOL-"- 

27X"  KM  GM  CM-2  KM  KM"^> 

FORMRT  ( SH  W ' 1 -4  • = 4 ■.  E 1 4 . 3 '•  ' 22X  > E 1 4 . 3/  > 

FORMRTC '/"CHOICE  OF  TWO  PRTHS  FOR  THIS  CRSE  -SHORTEST  PRTH  TRKE 
I N."  "FOR  LONGER  PRTH  SET  LEN=1.";' 

FORMRT'  / ' IX’  "LL'  - SX’  "LEVEL".  2X»  "RLTITUDE  TEMP",6X>  "PRES"'' 

1 7! . WH2O " ) 7X  ’ " NO  J " y OX " WGRS " / 9 
FORMRT'::/ "DISK  I/O  ERROR!.  I ERR  =".I3) 

END 


NO  ERPURS**  PR'OGRRM  = OlOSb  COMMON  = 055:3:3 


FhioE  yti01 


FTN4  COMPILER:  HP24ir? 


SEPT.  1974) 


FTO.  L 

PFOGRRM  SECTS  'S-'nT) 

D I MENS  I ON  TRHN 1 4 1 ■ , 1 Iii  B ' 1 44  ' » I S I ZE  < J:  ■ . NmM  1 ' 5 
D 1 MENS  I ON  NAME ' “T  < , T fill 21"  WH20 (.  40  ' !•  WO 3 ' 40  ' > WGMS  ( 40  > 

O'  iMMON  TEMP  (40 - PRES  ^ 40 ) < V 1 j V2 » DV » I RhU  . LMHX . RHZ2  (20  ' 
t.  OMMON  RANGE' UW^ 4 0.0  ' 

r OMMON  Z ( 34 P ' 34  ' . T 34  < EH  ( S j 34 ) ? UH  ■ 34 ’ M - NL  ' PE  < CM  > L 3 ? P I 
COMMON  WO  ( 34  " HZ  1 ' 34  ■ ■.  HZ2  •:  6 ) » TM 1 0 > . VH ' S ' • N (,  3 " E •:  A ) 

C OMMON  C4  ■:  1 33  » C S ' 1 S " HHHZE  ■'  34  ' < VK >'  45  "HZ'  45  ' 

COMMON  C7H'.45  :'  > I VP':  4u  > « ALT  ' 40> ' IPRM 
C OMMON  X 1 » X2 . H 1 ' H2 . N > NP  ? HP  1 ? NP2 » TX 1 » TX2  j (11 , VN 1 > VN2 
C OMMON  L • JP » L 1 ' L2 ' PS ' TS , X > REF  j I P ' LBR ' EV ' I MP , K2 1 H j DS 
C OMMON  PS  I , SPH I W HET ' T HET A » PH  I j BET  j BET  A , SALP ' SR 
COMMON  HM I N ' LL » HM ' RX 1 CA , RN  j LEN  > SUMA , AB ' AL AM 
C OMMON  I XV . V I S » I V 1 , I V2 ' I HAZE  > I TYPE ' ANGLE  - F I LL ' 1 392  > 

COMMON  T 1 . T2  * TM I N ' P 1 , P2 » PM  I N 

EQU I VALENCE  ' T AU ' AHZ2  ' ' ( WH20  j WW  ( 1 ' 1 .> ' W03  ? WW  a j 2 ) ) 

E OU I VALENCE  ' WG A S ' WW ' 1 ? 3 ) ) 

DATA  NAME  / 2HSE  j 2HGT  ? 2H5  -•  4 NAM1/2HTR^  2HAN,  2H1 
^ BEGINNING  OF  TRANSMITTANCE  CALCULATIONS 
ISIZE' !■  = 150 

I F ■ I RAD . NE . 0 .)  CALL  GREAT ' I DCB  j I ERR ' NAM  I , I S I ZE » 3 ) 

IF'  lERP.EQ.-l  :•  GO  TO  110 

IDV  = DV 

IV  = I VI -IDV 

LMAX  = LBR 

LOOP  = 1 

INDEX  = 0 

IF'  IRAD. EQ.  1 ' LOOP  = LBR 

IV  = IV+IDV 

INDEX=INDEX+1 

IF' INDEX. NE. 1'  GO  TO  52 

IF^ JP.NE.0)  GO  TO  52 

WRITE(b»422) 

CONT I NUE 
DO  95  LL=1»L00P 
DO  53  K=1.8 
WOO  = VH(K) 

IFcLL.GT. 1)  W<K>  = W' K >-WW(LL-l 5 K> 

TX(K;'  =0.0 

CONTINUE 

TX(l)  = 1.0 

ICOUNT=ICOUNT+1 

SUM=0. 0 

V=IV 

I=(IV-350>/5+l 
IF(IV.LT.670 ' GO  TO  72 
IF'  IV.LE.3000>  GO  TO  bl 

^ MOLECULAR  oLHfTERING 


Cb=9. 807E-20*(V#  + 4. 01 1 

TX(3)=C6*W<3> 

SUM=SUM+TX<3) 
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yy57 

IF  ( IV. LT. 13000  ' GO  TO  72 

yyss 

c 

‘3059 

r>  ♦ * + ♦ 

WRTER  VAPOR  CONTINUUM  18  MICRON  REGION 

00t'0 

c 

006 1 

ti 

IF'  IV.uT. 1350  ' GO  TO  62 

0062 

T V ( 2 ^ ' 4 . 1 8+5578 . O+EXP  ‘ -7 . 87E-  3#  VV)  *N  ( 2 ) 

0063 

GO  ro  t>6 

0064 

K J 

IF'- IV.  LT.  2350.:'  GU  TO  68 

0065 

c 

0066 

C » 

■ WATER  VAPOR  CONTINUUM  4 MICRON  REGION 

0067 

c 

0068 

X I = < (■'  - 2 3 5 0 . 0 .'  ■ 50.0+1.0 

0069 

HO  63  NH=1>15 

0070 

xh=xi-float(nh;' 

0071 

TX(  2 .'=C5(NH> 

010172 

IFO<H‘  64  >65 -63 

0073 

CONTINUE 

0074 

64 

TX ( 2 > =TX ( 2 ) +XH* ( C5 ( NH> - C5 ( NH- 1 > ) 

0075 

65 

TX(2)  = TX'.:2)+W(8;' 

0076 

66 

SUM=SUM+TX(2) 

0077 

IFdV.LE.  1350.OR.  IV.GT.2740)  GO  TO  72 

0078 

c 

0079 

NITROGEN  CONTINUUM 

010801 

c 

0081 

68 

IF  'IV.LT.20S0)  GO  TO  72 

0082 

K 4 = I - 3 4 6 

0083 

TX‘  1 )=C4<K4>+l'Kl  ■ 

0084 

SUM=SUM  + TX(1  :' 

0085 

c 

0086 

AEROSOL  EXTINCTION 

0087 

c 

0086 

72 

ALAM=1 . 0E+4.-V 

0089 

XX=0.0 

0090 

'I'V  = 0,0 

0091 

IF  (IHA2E.EQ.0. ) GO  TO  90 

0092 

DO  38  N = 1 » 4 4 

0093 

XD=ALAM-VX(N ' 

0094 

IF(XD>89>88»88 

0095 

'6S 

CONTINUE 

0096 

y 9 

XX=  ( C7  ( N )-C7(  N- 1 ;•  ' #XD/  ( VX  ( N > -VX ( N- 1 ) > +C7(  N ) 

0097 

YV=(C7A(N)-C7A(N-1 ) )#XD/(VX(H)-VX(N-1 > )+C7A(N 

0098 

90 

TX(8>  = VY*W(4:' 

0099 

TX(4>=XX*W(4:.' 

0100 

SUM=SIJM  + TX(4> 

0101 

TX(7>=SUM 

0102 

DO  94  K=l»8 

0103 

IF  (TX(K).EQ.0.0:'  GO  TO  92 

0104 

IF  CTX(K>.GT.20. ) GO  TO  93 

0105 

tx(;k:'=*exp<-tx(K);' 

0106 

GO  TO  94 

0107 

92 

TX(K>*1 . 0 

0108 

GO  TO  94 

0109 

93 

TX(K)=0. 

0110 

94 

CONTINUE 

0111 

TX(8)=*1.0-TX(8> 

0112 

AB=*1.-TX<7) 

1974) 
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eu3 

0114 

0115 
01  IE 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 


IF  IV.  EQ.  IVl . OR.  1 V.  E".  IV2)  flB=»0. 

SUMR*SUMfl+PB+DV 

IF<JP.EQ.0.RND.LL.EQ.  1 ' WF  ITEC6. 423)  IVp  RLRf'MX(7)  - . 

1F=1.4)»TX<8>pSUMH 
I F ( I RRD . NE . 0 > TRRN 1 (L L = TX •:  7 ) 

95  CONTINUE 

Thu (INDEX)  = TX'7' 

TRRNlCLMRX+l)  « IV 
II  * <Lf1fiX+l..'*2 

iF(I  RRD . NE . 0 ■ CRI,  1.  WR I TF ' I DC8 » I ERF » TFRN 1 p I L ) 

IF<IEPR.EQ.-1>  GO  TO  110 
IFaV.GE.IV2'  GO  TO  96 
GO  TO  50 

96  NRITE(6p496) 

DO  97  LL=1pLMRX 
F R C = W W (.  L L p 3 ) 

WHZ'O'LL)  = WU(LLp5)/FRC 
WOSCLL)  = WW'LLpG>/FRC 
NGRS'.'LL..'  = FRC 
C 

TEMPORRRY  PR  INI  OUT 
C 

97  UP  I TE  ( 6 p 498  > LL  p LYR  ( LL.  ' p RLT  ( LL  > p TEMP  ' LL ) p PRES LL ) p WH20 LL ) 
IpUOS'^ldpWgrscll:' 

I PPM  = 1 
NRME<3'  = 2H6 

if<;irrd.eq.0j  go  to  102 

IF'  WRITF<IDCE:p  IERR»  IEUFp  -1  )>  1 10p  101  p 101 

101  IF'XLOSEUDCBp  lERR))  110p102p102 

102  CRl.L  EXEC <8.  NAME) 

CRLL  MRDR 

110  WPITE(6p500>  IERR 
STOP 

422  FORMAT  •::/.'///■■  FREQ  WRVELENGTH  TOTRL  " 

I"  N2  CONT  H20  CONT  MOL  SORT  REROSOL  REROSOL  INTEGRRTED" 
2/1Xp14HCM-1  MICRONSp4<4X5HTRRNS>p4Xp "TRRNS  RBS  RBSORP" 
3”TI0N” . ) 

423  KORMRT  '::I6p7F9.4pF10.2> 

496  FORMRT' ///IXp 'LL"p3Xp "LEVEL" p2Xp "RLTITUDE" p 3Xp "TEMP"p6Xp "PRES"p 
17Xp "NH20"p7Xp "W03"p8Xp "NGRS"/) 

498  FORMAT < I3p I6p  3F10. 2p  2Xp  3E1 1 . 3) 

500  FORMRT</"DISK  I/O  ERRORp  IERR  ="pI3) 

END 


♦ # NO  ERRORS^"-  PROGPRM  = 01602 


COMMON  = 05533 
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yyei 

FTN.l 

PROGRRM  SEGT6  (5’90,'i 

0003 

D I MENS  ION  I DC  b a 44  ■ , JDC  B ' 1 44  ' » I S I ZE ( 2 ' • NRM 1 k 3 ' » NRM2 ( 3 ) 

0004 

D I MENS  I ON  NUN ( 6 > ? KPTS ( 3 » 40  ' » CON ( 6 > » SPEC  < 6 » 2 ) > Z < 2 > 

0005 

DIMENSION  NfiNE(3.MPP':9.MTT(9>»TRRNl(4r'jTRRN2<41> 

0006 

DIMENSION  eUF(63) 

0007 

COMMON  TEMP ( 40  j , PRES ‘ 40  > , V 1 > V2 » D V > I RRD » LMRX  > TRU ( 2 1 > 

0008 

COMMON  WH20  ( 40  '■ » N03  (40  ' !>  WGRS ' 40 ) < VV  ( 250 > ■«  RK  ( 9 ? 250  > 

0009 

INTEGER  SPEC 

0010 

DRTR  CON  -'0.0>  .330E-3i0.0j75E-7»  1.6E-6. 

0011 

DRTR  NRME  • 2HSE»  2HGT’ 2H7  , NRM  1 ■ ^HTR i 2HRN 2H 1 ■■ 

0012 

DRTR  NRM2  '2HTR,  2HRN.  2H2 

0013 

C 

0014 

c*-*- 

• PERD  CONTINUUM  TRU' 3 FOR  PRD  I RT I ON  CRLC 

0015 

c 

0016 

lOPTN  = 0 

0017 

IF'  IPRD. EQ. 0 ' GO  TO  70 

0018 

I F ( OPEN I DCB  < I ERR  ^ NRM  1 - I OPTN ''  ' 500  j 60  ^ 60 

00 1 9 

60 

: F '•  RERDF  ( I DCB . I ERR . TRRN 1 > > 5O0 ,61  >61 

0320 

61 

I F ' RERDF ( I DCB  > I ERR . TRRN2  > > 50O  ? 62 » 62 

0021 

62 

RVl  = TRRNKLMRX+n 

0022 

RV2  = TRRN2CLMRX+1  ■ 

0023 

70 

CONTINUE 

0024 

IB  = 1 

0025 

ISIZEa;  = 50 

0026 

I F ' CRERT ( JDCB » I ERR ’ NRM2  * I S I ZE » 3 > ) 50O  > 95  > 95 

0027 

C 

0028 

C***^-< 

■ INSURE  THRT  (V1<V2)  RRE  HI  THIN  TRPES  RRNGE 

0029 

c 

0030 

95 

RERD(5.9O0)  DVM 

0031 

I F ( DVM . LT . 0 . 005 > DVM=0 . 005 

0032 

WRITE(6»938>  DVM 

0033 

nSPEC  = 6 

0034 

C 

0035 

PERU  TRPE  BLOCK  INTO  DISK  FILE 

0036 

c 

0037 

REWIND  8 

0038 

RERD ( 8 » 900 ) VM I N < VMfiX :•  NPT 

0039 

IF' V2.GT. VMIN;  go  TO  105 

0040 

100 

WR I TE ( 6 > 942 ) V 1 » V2 , VM I H j VNRX 

0041 

STOP 

0042 

105 

IF(V1.GE.VMRX>  GO  TO  100 

0043 

120 

IFC'/l.GE.VMIH)  GO  TO  122 

0044 

WRITE (6. 930)  Vl^VNIN 

0045 

VI  =«  VMIN 

0046 

122 

IF(V2.LE. VMRX>  GO  TO  123 

0047 

WRITE(6i931>  V2»VNRX 

0043 

V2  = VMRX 

0049 

123 

CONTINUE 

0090 

C 

0051 

C*-^***  PERD  PiT  VRLUES  FROM  UISK  FILE 

0052 

C 

0053 

RERDC  8 j 90 1 ) ( PP  < K ' » K>  1 j NPT ) 

0094 

RERD(8»901)<TT(K>  jK«=1»NPT) 

0055 

c 

0056 

DETERMINE  IHTERPOLRTION  POINTS  FOR  ERCH  LRYER 
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0057 

C 

0058 

CHLL  PTPTS(PP-  TT-LMRI  SKPTS) 

0059 

C 

0060 

C * + 

• RERD  IN  WAVENUMBER  BLOCKS 

0061 

c 

0062 

VCHKl  = VI -10. 

0063 

VCHK2  = V2+10. 

0064 

ILP  = 1 

0065 

125 

NUMl  = 1 

006fc 

R E R D ( 8 ' 9 0 0 ) V R j V B 

0067 

30  128  M=1.MSPEC 

0068 

P E R D ' . 8 ? 9 0 2 . ' 8 P E L ' M ' 1 ' » 8 P E L v M i>  2 ) > N U M 2 

0069 

NUMCM)  = NUMl 

0070 

NUMl  = NUM1+NLIM2 

0071 

NMIN  = NUM(M' 

0072 

NMRX  = NUMl-1 

0073 

I F <.  NMRX . GT  . 25u WP  1 TE  « 988 ) VR  j NMRX » M 

0074 

riQ  127  N=NMIN,NMRX 

0075 

127 

READ  ( 8 - 903 ) VV  ( N :■  ? ' RK ' f'  - N < K=  1 ? NPT 

0078 

128 

CONTINUE 

0077 

IF' VR.GT.VCHKl)  GO  TO  129 

0078 

WR I TE  6 » 933 ) VR , HUM (N  • , N=1  > 6) » NMRX 

0079 

GO  TO  125 

0080 

129 

I F ' V R . b E . V C H 1 2 > G 0 TO  16 0 

6081 

L 

0082 

C^^**4*  CfiLCULfiTE  TPfiHSMISSION 

0083 

C 

0084 

WR I TE ( 6 , 906 ) NUM ( 1 ' * NUM ( 6 > 

0085 

IFdLP.EQ.  1)  WRITE'  6-932) 

0086 

ILP  = ILP+1 

0087 

V = VV(NMIN) 

0088 

VO  = V 

0089 

N = 0 

0090 

135 

IF(V.GE.Vl)  GO  TO  136 

0091 

N = N+1 

0092 

V = V0  + FLORT(N'*DVM 

0093 

GO  TO  135 

0094 

136 

N = 0 

0095 

VO  = V 

0096 

140 

N = N+1 

6097 

RRD  =0.0 

0098 

RRDl  = 1.0 

0099 

FRCl  = 0.0 

0100 

C4***4  CONTINUUM  THU’S  FOR  RHDIHNCE  CHLC 

0101 

IF' IRRD.EQ.O  ■ GO  TO  1404 

0102 

1 401 

IF(RV2.GE.V>  GO  TO  1404 

0103 

DO  1402  LL*1-LMRX 

0104 

1402 

TRRNKLL)  = TRflN2(LL> 

0105 

RVl  = RV2 

0106 

IF'  RERDFdDCB,  lERR- TRRN2  :- ) 500-  1403-  140 

6107 

1403 

PV2  = TRflN2(LMflX+l  ■ 

0108 

GO  TO  1401 

0109 

1404 

CONTINUE 

0110 

DO  150  LL=1-LMRX 

0111 

DIST  = WGflS(LL) 

0112 

CONd)  * WH20CLL) 

i.  UUrU  1 ^ * WMiLKLLJ 

I 

■ 
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0113 

COHO'  = UOS'OL) 

0114 

FhC2  = 0.0 

0115 

PBhR  = PRESCLL' 

0 1 1 t' 

IF' PBflR.LT.75.  ' PBflR=75, 

0117 

no  142  M=1,NSPEC 

0110 

IF' flP'  1 > NUMc:M>  :' . EO.  0. 0)  GO  TO  142 

0119 

,H  = NUM'iM' 

0120 

1405 

'v'Vl  = VV-iMl' 

0121 

VV2  = VV(I11  + 1 ■ 

0122 

IFO'.LE.VV2)  GO  TO  1406 

0123 

Ml  = Ml  + 1 

0124 

GO  TO  1405 

0125 

1406 

CONTINUE 

0126 

141 

no  1411  1=1,2 

i 

0127 

N1  = Ml+I-1 

0128 

Y0  = flK<KPTS':2,LL ',N1  ' 

0129 

FT  = FKY0.flK' KPTSa'LL>,Nl),TT';.KPTS'::2<LL). 

, TTCKPTS':.  1 , LL/ ) 

0130 

l-TEMP'.LD) 

0131 

FP  = F1(Y0,RF(KPTS':3,LL),N1),PP<KPTS':2,LL)' 

',PP(KPTS^3,LL:'.:' 

0132 

1 ,PBflP ) 

0133 

flkP  = FT+FP-Y0 

0134 

IF(flKK.LT.0>  flKK=0 

1 

0135 

IFCPBflP.LT.  75.  l)flKK=flKl:>TEMP(LL)/216. 6 

! 

0136 

IF(W(N1).EQ.V>  GO  TO  1412 

0137 

Z'l)  = RKK 

0138 

1411 

CONTINUE 

0139 

hKK  = F 1 ( 2 (1 1 > , Z 'i  2 , V'v'  1 5 VV2 , V ) 

0140 

1412 

FflC2  = FflC2  + RKK*CON':;M)*niST 

0141 

142 

CONTINUE 

i 

0142 

FflC2  = FflC2*1.0E5 

0143 

IFURflli.EQ.O)  GO  TO  149 

0144 

TRN  =0.0 

0145 

IF':FftC2.LT.20.  :•  TRN  = EXP':-FflC2) 

0146 

PLANK  = BLflM'OEMP';:  J)^  V) 

0147 

PRD  = TRN*RRD+<1.0-TPH)*PLRNK 

0148 

FRC4  = FiaRflNl(LL)*TRRN2':LL),RVl,RV2i' V) 

0149 

FRC5  =1.0 

0150 

IF' LL.NE.LMRX)  FRC5  = F 1 TRRN 1 CLL+ 1 ) , TRRN2 

;ll+i),rvi,pv2,v) 

0151 

RRDl  = RRDl  + (FRC5-FflC4::'+PLRNK 

0152 

149 

FflCl  = FRC1+FRC2 

0153 

150 

CONTINUE 

0154 

PAD  = RRD+RRIH 

0155 

TRRN  =0.0 

0156 

IF(FflCl.LT.20. ) TPRN  = EXPC-FRCl) 

0157 

BUF(IB)  = V 

0138 

BIJF<IB+1>  = TRRN 

0159 

BUF<IB+2)  = PRD 

0160 

IB  = IB+3 

0161 

IL  = <IB-1)*2 

0162 

IFUB.NE.64;  GO  TO  155 

0163 

IB  = 1 

0164 

I F ( WR I TF ( jnCB , I ERR , BUF , I L ) ) 500,1 55 , 1 55 

0165 

155 

HP  I TEC 6, 934/  V, TPRN » PAD 

0166 

V = 70  + FLOAT <N)*DVM 

0167 

IF^V.GT.V2)  GO  TO  160 

0168 

IF^V.LE.VB)  GO  TO  140 

\ 

I 
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0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 


N = A-1 
GO  TO  125 

1 60  I P ' WR I TF  i.  JDC  E « I ERR » BUF  > I L ' ' 500 » 1 65  j 1 65 

165  IF' WRITF(JDCE- IFPR-BUFj -1' ) 500,170-170 

170  I F ' I RAD . EO . 1 IF'  PURGE  ■:  I DC  B , I ERR , NAN2  ;■ ,'  500 ,175,175 

175  IF CLOSE ' JDCB , I ERR ) ) 500 ,180,180 

180  CALL  EXEC C 8, NAME  ' 

CALL  MRDA 

500  UR  I TEC  6, 940)  I ERR 
STOP 

900  FORMAT' 2F 10. 2, 15) 

901  FORMAT' 9F1 0.2) 

902  FORMAT C2A2, 15) 

903  FORMAT' FI  2.2- 9E12. 6) 

906  FORMAT C " NH I N - NMAX  = " ,2110  ' 

930  FORMAT-:  ."VI  , F10. 2, 5X*  "TOO  SMALL,  RESET  T0",F10. 

931  F0RMAT'..//"V2  =" , FIO.  2, 5X,  "TOO  LARGE,  RESET  TO",F10. 

932  FORMAT' FREQ  TRAN  RAD",-) 

933  FORMAT' "BLOCK  SKIPPED,  V =",F5.1,"  UAVENUMBERS" , 716) 

934  FORMAT ' F 1 0 . 2 - F 1 0 . 4 , 2E 1 0 . 3 ) 

936  FORMATC2I10) 

938  FORMAT'. - "MEDIUM  RESOLUTION  DV  =",F5.1,"  UAVENUMBERS" 
940  FORMAT' /.•■■"DISK  10  ERROR,  I ERR  =",I3) 

942  FORMAT-’ "TAPE  OUT  OF  RANGE  OF  DATA"/"V1  =",F7.1,"-  V2 
1.F7.1-",  VMIN  =",F7.1,",  VMAX  =",F7.1) 

END 


**  NO  ERRORS)-*  PROGRAM  = 02397 


COMMON  = 05450 
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PRGE;  0001  FTH4  CUNFILER:  HP24177  ‘.SEPT. 


0 1 9":* 

FUNCTIUN  BLRM'  T>  V;' 

0200 

Bl  HM  = 1 . 1989E-lb*'./**5.  EXP' 

1 . 43;38^V.-  T- 1 . ' 

0261 

RETUPN 

0202 

END 

tJO 

ER'FmjRS  *•+  PPOGRFiM  = 0’06‘^'U 

L ONMON  = 06000 

PhGE  0001  FTM4  COMPILER:  HP24177  'CSEPT 


I 01S^  FUNCTION  FHYl  jY2.:;1jX2jK) 

j 01SE  FI  = i' 1 + ( Y2-V  1 ;■  * (X-X I (X2-X 1 

j •eiS?  RETURN 

I 0198  END 


1974) 


. 1974;' 
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ytiOl 

tnjCi 

0003 

0004 

0005 
000 1. 
0007 
0000 
000? 
0010 
001  1 
0012 

0013 

0014 

0015 
00  It 
0017 
001c 

0019 

0020 
0021 
0022 

0023 

0024 

0025 
002t 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
0o3t 

0037 

0038 

0039 

0040 

0041 

0042 
0S43 

0044 

0045 
0O4t 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 


PhGE  0001  F7N4  COMPILER:  HP24177  ('SEPT.  1974) 


FTH.  L 

P F' 0 ij R AM  3 E G T 7 ( 5 ^ 9 0 ) 

CUflMON  TEMP  (40 ) j PRES ' 40  ' ? V 1 , V2  ^ DV  > I RhU  . LMhX  ? TRU  ( 2 1 ' 
CijMMON  WH20  (40 ) . W03 ( 40  > . WGRS  (40 ) « VV  ( 250 ) , RK  ( 9 « 250  ' 

I' I MENS  I OH  V'.  1001  ■'  - NfiMS'/S  ' 

D I M E N S I 0 H J D C B (1  4 4 ' » B U F'  ( 1 3 ) j T R ' 2 ' 1 0 5 0 > 

EQU I valence  ' WH20 , TR ' 1 • 1 ) ) 

I'lATA  NAM2- 2H  FRj  2HANj  2H2  ' 

If  = 1 
NUMV  = 0 
lOPTN  = 0 

I P OPEN  JDCB ' 1 EPF: , NAM2  < I OPTN  ■ ) 1 00  - 1 0 - 1 0 

10  I F P E A D F J D U B ^ I E F'  R » B U F ' 1 2 8 j L E N ' > 1 0 0 ’ 1 1 j 1 1 

1 1 IF'  LEN.EQ.-l  ' GO  TO  20 
LEN  = LEN't 

NUMV  = NUMV+LEN 
riO  15  1 = 1,  LEN 
V'  IK)  = BUFa*3-2  :• 

1 , IK  > = BUFa*3-l 
TR'2'IF.  ) = BUF  a 
IF  = IF+1 
15  CONTINUE 

GO  TO  10 
C 

'^OLD  IN  CONTINUUM  RE8IJLTS 
C 

20  IF'.  RWNDF ' JDCB , I ERR ) 1 00 ,21,21 

21  [MAX  = I-l 

'JO  25  1 = 1,  NUMV 
CAL.L  CALC2(V'.  I > , T,' 

25  TR'1,I'  = T*TR(1,I) 

IF-:  IRRD.  EQ.  0.'  TR':2,I.'  = 0.0 
UP  I TE - t , 902 ) 

I 1=0 

Du  70  1=1, NUMV 
11=11+1 

WF:  I TE  ( t , 903 ) V ( I , TR  - 1 , I > , TR  2,1) 

BUF(II+3-2)  = VU;:- 
BUF(II+3-l>  = TRa,I) 

BUF'..  II+3)  = TR(2,I) 

IF-II.NE.21)  GO  TO  70 
IL  = 126 

I F ■ WR I TF  < JDCB , I ERR , BUF , I L ) ) 1 00 , 55 , 55 


c:c 

70 

11=0 

CONTINUE 
IL  = 11*6 

30 

100 

IF (WRITE 'JDCB, I ERR, 
STOP 

WRITE-:  6, 904)  I ERR 
GO  TO  80 

BUF, IL; 

')  100,30,80 

902 
90  J 

FORMAT'/'"  V 
FORMAT -:F6.  1,2F8.4  - 

TRANS 

RAD" /> 

904 

FORMRT- //"DISK  I.  O 
END 

ERROR, 

lERR  =", 13) 
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HHl'l  1 

yOO  V 
you  4 

yyyY 
yyyt 
yyyr 
yooc 
000'- 
0010 
001 1 
001c 

00  1 5 
0014 

001  j 

00  It 
0017 
00  IS 
00  1 9 
00  £0 
0021 
0022 

0023 

0024 

0025 
002t 
0027 
002S 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
003t 
0037 

0036 

0039 

0040 

0041 

0042 

0043 

0044 
. 0045 

0046 

0047 
0043 

■ 0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 


F hGE  0001 


FTN4  rriNPILER;  HP24177  -SEPT.  1974) 


F rN«L 

iUBPOUTIME  POim  ■ I'M- tMlP?  TX,  IF' ' 

C OMMOM  TEMP  ' 40  - PPES  '■  40  ■ < V 1 ^ V2 ' DV ' I PhD  - L I'Ih:  : < HHZ'2  ' 20  ' 

C OMMOM  PhMGE  ’ UW  ''40.6  ' 

L iJMf'D.iN  2 ' 34  ' ' P' ' < T ' 34  ' ? EH ' 3 « J4  ’’ « UFl ' J4  ' M < ML  < PE  » t W < L l.i  ^ P 1 
DIMENSION  TX' 10> 

F * *f*--**^**-k->r*^^**  * ii:  i *r  * ■>  t ^ iHi  *■  * * ■* 

C VMEPOUTINE  POINT  COMMJTES  THE  MERN  PETPhUVE  INDEX  hBuVE  MND  BELON 

C R GIVEN  rltituhe  RNL  INTEPPOLRTES  ET.'PONENTIRLLV  to  detepmine  the 

C EOLUVRLENT  RBSOPBER  RHOUNTS  RT  THRT  RLTITUDE. 

r 

C I'3  THE  HEIGHT  IN  QUESTION 

C T:,',7..'  RND  'i'N  RPE  THE  NEhN  PEFPRlTIVE  INDICES  RBOVE  RND  BELON  X 

C N I:  the  LE'.'EL  integer  CORRESPONDING  TO  X OR  THE  LEVEL  BELON  X 

C NF  =1  IF  X COINCIDES  NITH  MODEL  RTMOSPHERE  LEVEL  .IF  NOT  NP  = 0 

C T':'  l-6>  RPE  RBSOPBER  RMOUNTS  PER  fTM  RT  HEIGHT  X 

N=NL 
NP  = 0 

I:-  'X.LT.O.O'  X^O, 

I H * !*. . GT  . Z ^ NL  '* GU  T *J  4 

DO  1 I=1,NL 

N=I 

IF  'X-Z'  I.'-'  2.4,1 

1 CONTINUE 

2 J2  = N 

I = N-1 

FRC  = X-Z  ^ N .•  ( Z ('  -L:  ' -Z ' N ) 

F'X  1 = P < N ■'  4 ( p ( J2  > / P < N > > +*FRC 
T. ' 1 = T ' N ' ♦ ' T ' J2  > ^ T N ) ' +-.FRC 
T •'  9 = TXl 

t:  •.  10:'=  pxi 

N: . 1 =WH ' N > ♦ <: NH C J2 > ■ NH ■ N ' **FRC 

TX'  3)  = CO*F'X  1 TX 1 -4 . 56E-6*NX  1 + TX 1 ^CN 

T!'’'2-'  = CO*P'  J2  ' T 'I  J2  ■' “4 . 56E“6*NH '■  J2  T J2  ■>  *CN 

i:  :<  \ > = CO*P ■:  N ? / T •: N ) -4 . 56E-6*NH < N > + T N ' #CN 

TXF7>  = 0.5E-6-^'.:TX(2 '+TX(3)> 

■'N  = 0.5E-6*' TX':  1 >+TX'::3) ) 
iF  MP.EQ.0:)  GO  TO  9 
DO  3 L=1.7 
F =L 

IFa.EQ.7>  K=8 

I y EH ',  K . N > . EQ . 0 . 0 ')  GO  TO  3 

IF  '■..EH':  K.  N ' . GT.  1 000. 0;'  GO  TO  3 

TX ■ K ) =EH < K , N ' * C EH ■ K . J2  > /EH< K . N > > **FRC 

3 CONTINUE 
GO  TO  9 

4 NP=1 

IF  ':  I P . E C! . 0 ) G 0 r u 6 
DO  5 f;  = i.s 

5 TX' K :'-EH(K.N::' 

6 TX' 7 '=EH(?,N '-1. 

■■i  N=0. 0 

IF  ''N.GT.D  VN=EHC7.N-1.J-1.0 
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F'hGE  fiCi02  POINT  FTN4  COnP I LER : HP24177  'GEPT.  1974) 


0057  9 

■TONT I NUE 

O05t: 

IF  UP.EQ.  1 )NRITE'  6. 400  :■ 

' ! j N !>  NP ) TX  c 7 ' 

0059 

TX  - 7.)  = TXc7)  + l . 

0060 

)'H='r'H+l  . 

0061 

RETURN 

0062  C 

0063  400 

FORMR'I  FROM  RuINTr 

HEIGHT=",Flo. 

0064 

1.  INDEX  RBOVE  & BELON  X- 

'-2E11.4» ".IP 

0065 

2RM0UNT9  PER  M1  RT  X= " , 4E 1 0 . 3 • 

0066 

END 

* * N 0 E R R L 

iR:.*-^  PROGRRM  = 0076b 

COMMON  = 

t 
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> I 
RB 


CO  ro 


FRIjE  0001 


KTH4  COflf-'ILERt  HP2-177  ' SEPT.  1974) 


00£il 

FTN,  L 

0tik?2 

SUBROUTINE  PIPTS'  T'P,  I r<  IMhX 

tiOOj 

DIMENSION  PP':;9,.'  - IT  '.:9  :•  < f PTS'.: 

Oti04 

i;  OMMON  TEMP  ■:  40  ' , PRES ''  40  ' 

000? 

r 

OOOt 

C**  • 

PPOGRHM  WRITTEN  FOR  9 F<T 

0f«O7 

c 

00OS 

DO  r-..0  ,1=1. 1 NR 

0009 

PO  = PPES'..J) 

00  1 0 

TO  = TEMP(J,:' 

001  1 

c 

IF'.  I F RLC2 . OT  . 0 ) GO  T 0 50 

0 0 1 2 

r 

.^ChLC2  = 1 

0012 

l"*’ 

0014 

r ♦ 4-  + 

--  FIRST  OREL  RT  GIVEN  F , T — - 

0015 

»“ 

00  1 t' 

IF  ' PO.  GT,  PP'::5.-' . RND.  TO.  GT.  TT 

00 1 7 

I F ' P 0 . ij  T . P P .J , ' b 1 .1  T LI  5 

00  1 8 

FI  = 1 

00  1 9 

f 2 = 2 

0020 

f:s  = 3 

0021 

IFCT0. LE. TT(2)  > GO  TO  50 

0022 

11=4 

0023 

f:2  = 3 

0024 

K5  = 2 

0025 

CO  TO  50 

0026 

c 

' F ' . P 0 . G T . P P . 5 , ' ) G 0 T 0 1. 0 

0027 

K1  = 3 

0028 

NS  = 4 

0029 

F3  = 5 

0030 

IF^  TO.  LE.  TT':  5.' ) GO  TO  50 

0031 

K 1 =6 

0032 

K2  = 5 

0032 

K2  = 4 

0034 

GO  TO  50 

0035 

10 

K1  = 6, 

0036 

K2  = 5 

0037 

F . ■:  = 7 

00  38 

PM  ID  = . 54'::PP'::5)  + PP';  7)  > 

0039 

IF':P0.LT.PMID>  GO  TO  50 

0040 

11=8 

0041 

F2  = 7 

0042 

K3  = 5 

0043 

GO  TO  50 

0044 

15 

I F ' F' 0 . G T . P P ' 7 ) ) C 0 TO  25 

0045 

K1  =9 

0046 

K2  = 8 

0047 

K3  = 6 

0048 

I F ':  T0 . GT . TT  (.  6 ) GO  TO  50 

0049 

R 6 - ( T 0 ~ T T ( 6 ) ) * * 2 + 'I  P 0 ~ P P ' . 6 

0050 

R7  = (TO-TTL? ')*,r2+<P0-PP'::7 

0051 

IF'R6.GT.R7)  GO  TO  20 

0052 

FI  = 5 

0053 

K2  = 6 

0054 

F3  = 8 

0055 

GO  TO  50 

0056 

20 

IF(T0.GE.TT<8))  GO  TO  30 

LOChTE  I N 1 E RPULRT I ON  Pu I NTS 


^ 1 .>■  ifl  i>  1 I 

• ^ , 


1 i. 
-i  i--R  -.L 
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FTH4  CUNPILEP:  HP24177 


';SEPT.  1974) 


PRGE  0002  PTPTS 


0057 

R1  = S 

0050 

K'2  = 7 

0059 

K3  = 5 

0060 

GO  TO  50 

0061 

•-.c 

TMID  = .5 

*(T1  '■  7>t  TT(U)  :■ 

0062 

K.  1 = y 

0063 

K2  = 7 

0064 

K3  = 5 

0065 

IF' TO. LE. 

TMIIt)  GO  TO  50 

0066 

n =7 

0067 

K2  = S 

0060 

K3  = b 

0069 

GO  TO  50 

0070 

30 

K1  = 9 

0071 

K2  = .9 

0072 

K3  = 6 

0073 

50 

CONTINUE 

0074 

KPTS<  1 ' .1,' 

= K1 

0075 

KPTS';2'  J> 

= f2 

0076 

KPTS<  3>  J ..' 

= K3 

0077 

60 

CONTINUE 

0076 

RETURN 

0079 

END 

**  NO  ERRORS**  PROGRhN  = 0042b  COMMON  = OOlbO 


'.1- 


'‘-vyi 


y 
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f-'HGE 


I ■ ■ j ' ■ • . 

{ • . 

bi— — 

riH4  COMPILER:  HP24ir?  tPEPT.  1974) 


0001 

FRUL 

0002 

SUBROUTINE  IHLC 

2 ' T2 

0003 

COMMON  TEMP' 40.:' 

- - RES ''40  , V 1 « V2  j DV  < 1 RhD  < LMRX , TRU ( 2 1 ) 

0004 

VI  = VI 

0005 

VF  = Vl+DV 

0006 

DO  10  1=1 » 19 

0007 

I F ':  V . GE . V 1 . RND . 

V.LT.VF'  GO  TO  20 

0008 

VI  = VI  + DV 

0009 

10 

VF  = VF  + DV 

0010 

20 

T2  = TRU'  I ■ + ' 

F TRU':.  I + 1 > -TRIJF  I > ) ■-■DV + ' V-V I ;' 

001  1 

RETURN 

00 1 2 

END 

ChOE  ny01  FTN4  COMP  I L EP : HP24177  ‘ OEPT.  1974) 


0001  FTNi.L 

0 002  F U H C T 1 0 N H 2 1 H fl  ' 

0003  TEST  = 1.0  - FIBS' H ' 

0004  hPG  = H -SQPT'  1 . -HrT*2  ■ 

0005  RSIN  = hTRH'.RRG' 

OOOb  IF' TEST.  LT.  1 . E-y  RSIN  = 2 . 0*RTRN i:' 1 . O > 

0007  RETURN 

0009  END 


I RGE  0001 


FTN4  COMPILE?:  HP24177  'SEPT.  1974) 


0001  FTN>L 

0 0 0 2 F U N C T 1 0 N R C 0 S ' . R ' 

000  3 RPG  = SQRT(  1 . -R*  + 2 :'/R 

0004  RlOS  = RTRN(RPG) 

0005  IF' fi. LT. 1 . E-S)  RCOS  = 2. 0*RTRN( 1 . O) 

000t  RETURN  , 

0007  END  ^ 


PRCE  0001  FTN4  COMPILER:  HP24177  'SEPT.  1974) 


0001  FTN.L 

0002  FUNCTION  F(:R' 

0003  F = EXP < 19. 9766-14. 9595*R-2. 43882*R*R>*R 

0004  RETURN 

0005  END 


■.‘I 
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APPENDIX  B 
MRDA  FLOW  CHART 


B-1 


SKCT  1 


SKC.T  2 


sec;  r ;< 


t 


I 

! 


H-4 


SEG'I’  (i 


H-7 


SKl'i  T / 


APPENDIX  C 


UST  OF  SYMBOLS 


SFOMKN'r  1 


Al^\M 

AVW 

BET 


J5F.TA 


CA 

CO 

t Vv 

C7 

C7A 

DP 

DV 

MMIN 


VVavflcn^;tli  (/jim) 

Input  zcmUi  :uik1''  (ilcHrt-i'S)  (c-ompar'  will  m U • t<  \l 

Avnrani’  wavelength  used  in  relraelivt  ind<  \ i pri  ^ m 

Anftle  subtended  at  the  eartli’.s  eenlt  r .is  path  ir  .M  rsi  .id  ai  ' nt 
levels  (cf  p.  in  Fc(.  (H)) 

Total  anfjle  sulitended  bv  path  at  earth's  eenter  (compai'  ii‘  1 q. 

(‘J) 

Conversion  factor  from  def>rees  to  radian.-' 

Waveleni^h  dependent  coefficient  used  in  relraetivi-  index  i .x]>ression 
Wavelength  dependent  coefficient  used  in  refractive  index  expression 
Extinction  coefficient  for  aerosol  models 
Aerosol  alisorption  coefficient 
Dew  point  temperature  (^C) 

Wavenumber  increment  at  which  transmittance  is  calculated 
Minimum  altitude  of  path  trajectory  (km) 


HI 

HJ 

HZl 

IIZ2 

I 

IIIAZF 

IM 


» IP 


IPRM 

ITYPF 

IXY 

J 

JP 

,11 

Xi 

L 


Initial  altitude  (km) 

Final  altitude  (km) 

Aerosol  number  density  (no.  cm  ^)  for  2;J  km  visual  range' 

_3 

Aerosol  number  density  (no.  cm  ) for  3 km  visual  range 

Running  integer  used  as  altitude  (level)  indicator  and  frequency 
indicator 

Aerosol  model  indicator 

Parameter  used  when  reading  in  a new  atmospheric  model  (.see 
Section  .'3.2.  1) 

Indicator  for  u.sing  sulirouUne  POINT  to  calculate  refractive  inck'x 
only  (IP  = 0)  or  equivalent  absorln'r  amounts  also  (IP  ^0). 

Flag  for  segment  communication 

Indicator  for  tyiie  of  atmospheric  path  (see  Section  .">.  I) 

Parameter  for  terminating  program  and  cycling  indicator 

Running  integer  lor  altitude  identification 

Print  option  parameter 

Level  indicator  for  altitude  Ml 

Level  indicator  for  altitude  H2 

Index  for  levels 


C-2 


LEN 

M 

ML 

MODEL 

NAME 

NL 

NLDAT 

PO) 

PI 

RANGE 

RE 

RH 

RI 

SPHl 

SR 

Ta> 

TRACK 

TMP 

TT 

VH(K» 

VLS 

VX 

VI 

V2 

WHO) 

WO(l) 

X 

XI 
X2 
Y 

Zfl) 


Parameter  used  for  defininj?  longer  of  two  paths  (see  Section  5.1) 
Integer  used  to  identify  required  model  atmosphere 
Number  of  levels  in  radiosonde  data  input  (MODEL  = 2) 

Integer  used  to  identify  required  model  atmosphere  (see  Section  5. 1) 

ASCII  Name  of  next  segment  to  be  called 

Number  of  levels  in  model  atmosphere 

Number  of  layers  in  model  atmosphere  data 

Pressure  (mb)  at  level  I 

3. 14ir)92f>54  that  is  (tt) 

Path  length  (km) 

Earth  radius  (km) 

Relative  Immidity  (%) 

The  product  of  the  sine  of  the  initial  zenith  angle  and  the  earth 
center  distance  to  starting  altitude 

Sine  of  the  local  zenith  angle  at  a given  level  (cf  sin  6 ) 

Slant  range  (km) 

Temperature  (°K)  at  level  1 

Background  radiation  calculation  temperature 

Ambient  temperature  (°C) 

Ratio  273. 15/ (TMP  + 273.  15) 

Integral  of  the  equivalent  absorber  amounts  from  HI  to  level  1 
Visual  range  (km)  at  sea  level 

Wavelength  at  which  aerosol  coefficients  are  read  in  (^m) 

Initial  frequency  for  transmittance  calculation,  cm  ^ 

FMnal  frequency  for  transmittance  calculation,  cm  ^ 

_3 

Water  vapor  density  at  level  I (gm  m ) 

_3 

Ozone  density  at  level  I (gm  m ) 

Input  height  to  POINT  subroutine 
Earth  center  distance  of  level  I 
Earth  center  distance  of  level  I + 1 
Input  zenith  angle  in  radians 
Altitude  at  level  I in  km 


SKC.MKNT  2 


AHAZK 

AH  7.2 

ALT(l.L) 

ANGLK 

CO 

CW 

D 

[)S 

DZ 

KH 

EV 


.Xc-rosol  lumilH  r density  for  MODEL  1 
Aerosol  imml)er  density  for  Mf)I)KI,  2 
Altitude  at  level  Z(I.YIl(LL)) 

Input  zenith  an>4^’  (deftrees)  (compare  witli  6 in  the  t(“.\t) 

Wavelength  dependent  ccxifficient  used  in  relractive  index  expression 

Wavelength  de])endent  coefficient  usi-fl  in  refractive  index  expression 

Water  va[)or  amount  ([)r.  cm/km)  at  level  1 

Path  length  from  level  l to  Level  1 + 1 

Heijtht  increment  from  level  I to  level  I t-  1 

Equivalent  absorber  amounts 

Integrated  absorl)er  amount  from  level  1 to  level  I +•  I 


HAZE 

111 

H2 

HZl 

HZ2 

IHAZE 

IP 

IPRM 

ITYPE 

JP 

L 

LHR 

LL 

LYR(I.L) 

M 


Aerosol  number  density  (no.  cm  '*) 

Initial  Altitude 
Final  Altitude 

Aerosol  number  density  (no.  cm  '*)  for  2.'1  km  visual  rangi' 
Aerosol  number  density  (no.  cm  for  .a  km  visual  range 
Aerosol  model  indicator 

Indicator  for  using  subroutine  POINT  to  calculate  refractive  index 
only  (IP  = 0)  or  equivalent  absorber  amounts  also  (IP  / 0). 

F’lag  for  segment  communication 

Indicator  for  type  of  atmospheric  path  (see  Section  .5.  If 
Print  option  parameter 

Frequency  indicator  for  ozone  transmittance  calculation 
I.ayer  counting  parameter  for  slant  path  trajectory 
Ixivel  index 

Ijcvel  number  of  Lth  Layer  in  the  atmospheric  path 
Integer  used  to  identify  model  atmosphere 


C-4 


X 


NAMK 

NL 

NI*1 

PHI 

I>S 

PSl 

PPVV 

PRES(LL) 

}>T 

REF 

RN 

RX 

SALP 

SPHl 

SR 

TEMP(LL) 

TO) 

THETA 

TS 

TX(K) 

TX(9) 

TX(IO) 

TXl 

VIS 

WH(I) 

ww 

W2 

X 

XI 
YN 

za> 


ASCII  of  next  sepnent  to  be  called 

Number  of  levels  in  model  atmosphere  data 

Value  of  NP  for  altitude  111 

An^le  of  arrival  at  H2 

Total  pressure!  in  atmospheres 

Angular  deviation  of  path  from  initial  direction 

Partial  pressure  H2O 

Pressure  al  level  l.L 

Product  of  total  pressure  (atm)  and  the  square  root  of  273/T(M,I) 
Refractive  index  of  air  at  level  1 

Ratio  of  refractive  indices  of  air  above  and  below  a given  level 
Ratio  of  earth  center  distances  between  adjacent  levels 
Sine  of  angle  of  arrival  at  adjacent  level 
Sine  of  the  local  zenith  angle  at  a given  level 
Slant  range  (km) 

Temperature  at  level  LL 

Temperature  (*^K)  at  level  I 

Zenith  angle  at  a given  level  (in  degrees) 

Ratio  of  standard  temperature  (273. 0°K)  to  temperature  at  let  -1  I 

Equivalent  absorber  amounts  per  km  at  a giv  en  altitude  obtained 
from  POINT;  also  transmittance  values  at  a given  wavelength  for 
each  absorber  type  (K  = 1,  8) 

Total  transmittance  at  frequency  V 

Absorption  due  to  aerosol  only  at  frequency  V 

Refractive  index  of  layer  above  initial  altitude  HI 

Visual  range  (km)  at  sea  level 

_3 

Water  vapor  density  at  level  I (gm  m ) 

Equivalent  absorber  amount  from  observer  to  level  L 
Water  vapor  density  for  atmospheric  model  M at  level  1 + 1 (gm  m 
Input  height  to  POINT  subroutine 
Earth  center  distance  of  level  I 

Refractive  index  of  layer  below  input  height  from  POINT  subroutine 
Altitude  at  level  I in  km 


SKC.MKNT  3 


I, 


AL 

Al.T(LL) 

AO 

HKT 

BKTA 

CA 

I)S 

KV 

H 

HM 

MMIN 

III 

M2 

1 

IP 

I PHM 

JP 

K 

K2 

!. 

MIH 

I.l. 

LMIN 

LYK(LL) 

LI 


Kquivalent  absorber  amount  i)cr  km  at  level  L 
Altitude  of  level  LYR(LL) 

('onstant  A defined  in  Kq.  (10)  of  LOWTRAN3  Manual 

Angle  subtended  at  the  eartli's  center  as  ))ath  traverses  adjacent 
levels 

Total  angle  subtended  by  path  at  earth's  center 
Conversion  factor  from  degrees  to  radians 
Patli  length  from  level  1 to  Level  1 + 1 
Integrated  absorber  amount  from  level  1 to  level  1 + 1 

Altitude  at  which  calculations  are  being  made 
Estimated  tangent  heif^it  (km) 

Minimum  altitude  of  path  trajectory  (km) 

Initial  altitude 
Final  altitude 

Running  integer  used  as  altitude  (level)  indicator  and  frequency 
indicator 

Indicator  for  using  subroutine  POINT  to  calculate  refractive  index 
only  (IP  = 0)  or  equivalent  absorber  amounts  also  (IP  / 0). 

Flag  for  segment  communication 

Print  oj)tion  |)arameter 

Absorljer  indicator 

Cycling  parameter  for  downward  looking  paths 

Frequency  indicator  for  transmittance  calculation 

Total  Number  of  levels  transvorsed  in  the  path 

liCvel  index 

Layer  number  of  UMIN 

Lth  level  in  path 

Frequency  identilier  for  transmittance 
calculation 

Frequency  identifier  for  transmittance 
calculation 
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I 


1 


M 

NAME 

NPl 

NP2 

P(I) 

PRES(LL) 

rei 

Range 

Rt: 

REF 

RN 

RX 

SALP 

SPHI 

SR 

TO) 

TEMP(LL) 

THET 

THETA 

TMP 

TX(K) 

TX(9) 

TX(IO) 

VH(K) 

W(K) 

WW(K,L) 

X 

X2 

Y 

YNl 

YN2 

Z(l) 


Integer  used  to  identify  required  model  atmosphere 

Used  to  contain  name  of  next  segment  called 

Value  of  NP  for  altitude  HI 

Value  of  NP  for  altitude  H2 

Pressure  (mb)  at  level  I 

Pressure  at  level  LL 

Angular  deviation  of  path  from  initial  direction 
Path  length  (km) 

Earth  radius  (km) 

Refractive  index  of  air  at  level  1 

Ratio  of  refractive  indices  of  air  above  and  below  a given  h;vel 
Ratio  of  earth  center  distances  between  adjacent  levels 
Sine  of  an^e  of  arrival  at  adjacent  level  (cf  sin  a ) 

Sine  of  the  local  zenith  angle  at  a given  level  (cf  sin  0 ) 

Slant  range  (km) 

Temperature  (°K)  at  level  I 

Temperature  at  level  LL 

Zenith  angle  at  a given  level  (in  radians) 

Zenitli  angle  at  a given  level  (in  degrees) 

Ambient  temperature  (°C) 

Equivalent  absorber  amounts  per  km  at  a given  altitude  obtained 
from  POINT;  also  transmittance  values  at  a given  wavelength  for 
each  absorber  type  (K  = 1,8) 

Total  transmittance  at  frequency  V 

Absorption  due  to  aerosol  only  at  frequency  V 

Integral  of  the  equivalent  absorber  amounts  from  HI  to  level  I 

Total  equivalent  absorber  amount  for  entire  path 

Equivalent  absorber  amount  from  observer  to  level  LL 

Input  height  to  POINT  subroutine 

Earth  center  distance  of  level  I + 1 

Input  zenith  angle  in  radians 

Refractive  index  of  layer  below  initial  altitude  HI 

Refractive  index  of  layer  below  final  altitude  H2 

Altitude  at  level  I in  km 
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SEGMENT  4 


Ai;r(LU 

B 

beta 

CA 

OV 

E(K) 

H 

HM 

1 

IPRM 

IVl 

IV2 

K2 

L 

LBR 

LEN 

LL 

LLNDN 

LMAP 

LSTORE 

LYR(LL) 

12 

NAME(3) 

P(L) 

PH 

PMIN 

PRES(LL) 


Altitude  of  layer  L 
Blaekbody  funetion 

Total  angle  subtended  by  path  at  earth's  center 
Conversion  factor  from  degrees  to  radians 
Wavelength  increment 

Equivalent  absorber  amounts  per  km  at  height  H 
Height 

Estimated  height 

Minimum  altitude  of  path  tragectory 
Frequency  index 
Segment  path  indicator 
Starting  frequency 
Last  frequency 

Cycling  parameter  for  downward  path 
Running  index  for  layers 

I^ayer  counting  parameters  for  slant  path  trajectories 

Parameter  used  for  defining  longer  of  two  paths 

Running  index  for  levels 

Value  of  LL  at  HMIN 

Counting  variable  for  long  path  storage 

Counting  variable  for  layer  index 

Altitude  of  Lth  layer  in  path 

Layer  number  of  H2 

Used  to  identify  next  segment  to  be  called 

Pressure  at  layer  L 

Angle  of  arrival  at  H2 

Pressure  at  HMIN 

Atmospheric  pressure  at  layer  L 
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PS 


PSl 

RN 

SPHI 

SR 

T(L) 

TEMP(LL) 

TMIN 

TS 

TXl 

VI 

VH(K) 

W<K» 

WW(LL,  K) 

XI 

X2 

YNl 


Total  pressure  in  atmospheres 

Angular  deviation  of  path  from  initial  direction 

Ratio  of  refractive  indicies  of  Air  above  and  below  a given  level 

Sine  of  local  zenith  angle  at  a given  level 

Slant  Range 

Temperature  at  layer  L 
Temperature  at  layer  LL 
Temperature  at  HMIN 

Ratio  of  standard  temperature  to  temperature  at  level  L 
Refractive  index  of  layer  above  initial  altitude  H 1 
Initial  frequency 

Integral  of  the  equivalent  absorber  amounts  from  HI  to  level  I 

Total  equivalent  absorber  amount  for  entire  path 

Equivalent  absorber  amounts  from  observer  to  level  LL 

Earth  center  distance  of  level  L 

Earth  center  distance  of  level  L + 1 

Refractive  index  of  layer  below  initial  altitude  HI 


SEGMENT  5 


AB 

ALAM- 

C4 

C5 

C6 

C7 

C7A 

DV 

IDV 

I HAZE 

INDEX 

IPRM 

IRAD 

IV 

IVI 

IV2 

IP 

K4 

LL 

LMAX 

LOOP 

NAME 

NH 


SUM 
SUMA 
TAU  (INDEX) 
IRANI  (LL) 


Absorption  at  frequency  V;  also  average  transmittance 
Wavelength  (/xm) 

Absorption  coefficient  for  nitrogen  (~4/xm) 

Absorotion  coefficient  for  water  vapor  continuum 

Extenction  coefficient  for  molecular  scattering 

Extenction  coefficient  for  aerosol  models 

Aerosol  absorption  coefficient 

Wavelength  increment 

Frequency  increment 

Aerosol  model  indicator 

Counting  variable  for  frequency 

Segment  path  indicator 

Radiation  calculation  flag 

Frequency  of  calculations 

Starting  frequency 

Last  frequency 

Print  option  parameter 

Frequency  indicator  for  nitrogen  continuum  transmission  calculation 
Running  index  for  level s 
Number  of  layers  in  the  path 

Number  of  layers  for  low  resolution  radiance  calculations 

Used  to  identity  next  segment  to  be  called 

Frequency  indicator  for  water  vapor  continuum  transmittance 
calculation 

Sum  of  optical  thicknesses  of  absorbers  4 through  8 
Accumulated  integrated  absorption 
Transmittance 
Transmittance  of  layer  LL 
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TX(K) 

TX(5) 

TX(6) 

VH(K) 

VX 

W(K) 

WGAS(LL) 

WHZO(LL) 

W03(LL) 

WW(LL,K) 

XD 

XH 

XI 

XX 

YY 


Equivalent  absorber  amounts  per  km  at  a given  altitude  obtained 
from  POINT;  also  transmittance  values  at  a given  wavelength 
for  each  absorber  (K  = 1,4) 

Total  transmittance  at  frequency  V 

Absorption  due  to  aerosol  at  frequency  V 

Integral  of  the  equivalent  absorber  amounts  from  HI  to  level  I 

Wavelength  of  aerosol  coefficients 

Total  equivalent  absorber  amount  for  entire  path 

Atmospheric  gas  density 

H2O  density 

Oj  density 

Equivalent  absorber  amounts  from  observer  to  level  LL 
Wavenumber  interpolation  parameter 

Wavenumber  interpolation  parameter  in  H„0  continuum 
calculation 

Wavenumber  interpolation  parameter 

Aerosol  extinction  coefficient 

Aerosol  absorotion  coefficient  of  frequency  V 
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SEGMENT  6 


AK(K,N) 

AKK 

BUF(63) 

CDN(6) 

DIST 

DVM 

FACl 

FAC2 

FAC4 

FACS 

FP 

FT 

ILP 

IRAD 

KPTS(3,LL) 

MSPEC 

Mi 

N 

NAME(3) 

NMAX 

NMIN 

NPT 

NVM(M) 

NI 

PBAR 


Extinction  coefficient  reed  from  tape  for  Kth  pressure- 
temperature  point  at  frequency  VV(N) 

Interpolated  extinction  coefficient 

Disk  write  buffer  containing  frequency  plus  spectra.' 
transmission  & radiance  results 

Species  concentrations 

Optical  depth  of  a species 

MRDA  frequency  interval 

Log  transmittance 

Summing  variable  for  transmittance 

Log  transmittance  for  radiation 

Log  transmittance  for  radiation 

Intermediate  result  in  interpolation  of  AK(K,N) 

Intermediate  result  in  interpolation  of  AK(K,N) 

Integer  variable  for  printing  heading 
Radiation  calculation  flag 

Elements  in  P-T  matrix  used  for  AK  interpolation 
Number  of  species 

Index  for  locating  extinction  coefficient  on  tape 
Frequency  index 

ASCII  name  cf  next  segment  to  be  called 

Upper  limit  on  index  for  input  of  W(N)  and  AK(K,N)  for  a particular 
species 

Lower  limit  on  index  for  input  of  W(N)  and  AK(K,N)  fer  a particular 
species 

Number  of  points  In  the  pressure  temperature  matrix 

Element  in  AK(K,N)  & W(N)  where  the  Mth  Species  information  beg 

Index  for  locating  extinction  coefficient 

Pressure  variable  (PBAH  > 75  mb) 


PLANK 

PP(K) 

PRES(LL) 

RAD 

RADI 

RVl 

RV2 

SPEC(M,2) 

TRAN 

TRANl(LL) 

TRAN2(LL) 

TT(K) 

VA 

VB 

VCHK 

VCHK2 

VMAX 

VMIN 

VO 

VI 

V2 

W(N) 

VVl 

W2 

WGAS(LL) 

WH20(LL) 

W03(LL) 

YO 

Z(I) 


Black  body  radiation  from  BEAM 
Pressure  array  read  from  tape 
Pressure  at  layer  LL 
Radiation  result 
Spectral  radiation 

Frequency  corresponding  to  information  in  IRANI  (LL) 

Frequency  corresponding  to  informtion  in  TRAN2(LL) 

ASCII  abreviation  of  Mth  specie 
Total  transmission 

Buffer  used  to  read  continuum  radiation  from  disk  for  layer  LL 
Buffer  used  to  read  continuum  radiation  from  disk  for  layer  LL 
Temperature  array  read  from  tape 
Initial  frequency  in  taoe  data  block 
Final  frequency  in  tape  data  block 

Used  to  compare  lower  frequency  of  tape  data  block  with 
calculation  frequency 

Used  to  compare  upper  frequency  of  tape  data  block  with  calculati 
frequency 

Max  frequency  contained  in  tape 
Minimum  frequency  contained  in  tape 
Initial  calculational  frequency 
Initial  calculational  frequency 
Final  calculational  frequency 
Frequency  array  read  from  tape 

Used  in  interpolating  tape  input  frequencies  to  calculation  frequen 

Used  in  interpolating  tape  input  frequencies  to  calculating  frequen 

Gas  concentration 

Water  vapor  concentration 

Ozone  concentration 

Used  in  interpolating  AK 

Altitude 
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SEdMKNrr  7 
HUE  (I) 

IK 

IRAD 

I.EN 

NUMV 

TR(1,  IK) 

rR(2,IK) 

V(IK) 


I 


Buffer  containing  frequency  plus  specular  radiation  and  transmission 
results 

Running  index  for  frequency 
Radiation  calculation  flag 

Nunil)or  of  frequencies  read  from  disk  in  one  read 
Total  length  of  V(IK)  array 

Combined  continuum  and  spectral  transmission 
Combined  continuum  and  spectral  radiation 
Frequency  corresponding  to  TR(I,IK) 


C-U 


’W05P  AH  LPL  (If<P  JT  ,0I»tc>UT,TAPE2,T  A«»Ei»  ,T  APES  = INPUT) 
? DCTOFER  7S  mTTO<vh  iOOTFTEO  FOR  MROA 


DIMENSION  W(7),  GHIMFOJO),  S(5000),  ALPHA(50J0),  EO°(ScnO) 
DIMENSION  M0L(50DD  ,VSTD?(4'13)  ,STOR{9,3L4),P(li)),r(11) 
DIMENSION  CS2  (9)  , OMPr, s ( 2]  i , 6)  , JCALC(6)  , SO  (533 0 ) , Al P 0 C50 0 0 ) 
DIMENSION  TI(250,12> , TTI (250) ,SPECIE(7) ,0ENS(6) 

-OSICAL  LOGIC 

DATA  SPEC1E/34H20, TMr02, 2H03f  3HN20, 2HCO, 3MCH4  ,2 HD 2/ 


2EHINO  2 
IEDF  = e 
DE=rH=t, 001 
>1=3.14159 

f 

lur7 

>L0HEP=1.0E-23 


i/BLDCK-lO. 

3REIT=.06 


i 


i. 


A t !'  i' 


^ im,  L 


REA0(9,76)  NPTPTS 

76  -0PM  AT  (12  ) 

?EA0(5,>7)  (P(I) ,T=1 ,NPT>TS) 

77  -OPMAT(>)  (EIJ.O)  ) 

?EA0(5,>7)  (T(I) ,T=1 ,NPT>TS) 

>RIMT  62, (P(I) ,T=1 ,MOTPT3) 

«2  -OPMAK*  PRESSURE  = *,'5(2X,F7.2)/10X,5(2X,F7.2)  ) 

>PIMT  84,  (T(I)  , T = 1,Notot5) 

84  -ORMAT(*  TEMP£PArt)PP=’»,5<2X,F7.2)/UX,5(2X,F-^.2)) 
rF(E0F(5 ) .NE.3)  STO® 

2EAD  81,  (W(M)  ,M  = l,-») 

81  -Oi?MAT(7E10.3) 

>PIMT  63 

83  -0PMAT(3X  ,*M  AT  ER*  , 6X  , »C  D 2*  , 6X , *0  70N  E* , 7 X ,*N  20  •,  7X , * 00  » ,8X , »CH4*  , 
irX,*02*,4X) 

>RIMT  81,  (M(M)  ,M  = 1,"') 

READ  85,  (/!,  Vi2,0</,POiJN0 

85  -ORMAT  (6F10. 3) 

•PINT  87,  Vl,'/V2,nv,'>0lJND 


D-2 


87  ¥1  ¥?  I7¥«^,P18.3,*  0<MNI>  , Fit . II 

C 4MTE(10,iaa)  ¥l,//?,*«PT>TS 

4H-ITE VI, vv?, Mayors 

120  "0«HAT «2 (F10.2) ,15) 

tP'(Vl.GE.  VV?)  STOP  ?1 
C ><RITE  <10, 13))  (0(1)  ,t  = 1,¥PTPTS) 

VPITE  (i*,  131)  (P(I)  ,T-rl,N»TPT3) 
r,  ^RITE(10,130)  (T  ( I ) , Tsi , »(PTPTS ) 

■<Wfe«»,131)  (T(I)  ,T3l,f<BTPTS) 

130  -0R8AT((5  (FIO.2) ) ) 

131  =3R'1AT  ((5  (Flf . 2) ) ) 

C 

r.  .030  6ACK  POINT op^o  MORE  LINES  FROM  THE  TAPE 

999  V2=V1*VBL0CK 

[ pfV2.6T. VV2)  V2=V79 
V90T  = Vl-B0UNn 
yf8P=V2tBOUHO 
VVTDP=VV2*B0UNn 

c 

C ^E  ARE  NOW  READY  TO  OEAO  TAPE. 

C . , 

t =1 

Itt=l  - 

VE3F=C 
REMIND  2 

1 READ(2)TMIN,TMAX,NI0EC, (<TI(II,J1), Jl=l,12) ,ITI(Il)  ,T l=l,NIREC) 
JT0=I0CHEC(2) 

TF ( JT0)4,2 

H *RINT  89,GNU(I) 

69  -ORNAT(*  PARITY  ERRO®  ENDOUNTEREO  AT*,F12.3) 

;0  TO  1 

2 [F(E0F(2) )5,7 

5 IE0F=IE0F^1 

VEOF=NEOF*1 

»riNT  91,IE^,TMIN,TMAK,NIREC  

91  -ORHAT(*  END  OF  FIL^  ENCOUNTERED*, 15, 2F1  2.  3,1  5) 

IFtfNEOF.GT.2)  STOP  ?» 

;o  TO  1 
? iEOf=€ 

rP(THAX.LT.VBOT)GO  Yn  l 
- m -9  Ksrl,NIRtC 

IF(TI  (K,1).LT.VB0T)  r-n  TO  9 
5NU(I)xTI (K,1) 

5(1) =TI(<,2) 

ALPHA  (1)*TI(K,3) 

-;DP(1)*TI  (K,i») 

MOL <I)«ITI(K) 

H=HOL  (1) 

5NIN*S(I) *OEMS (M) 

[F(SMIN,LE.l.9E-2‘.)  00  TO  9 
IF (GNU(I) .GT.VVTDP) GO  TO  11 
I=I»1 

9 CONTINUE 

C 

C MAXIHUM  NUNBER  OF  LINF^ 

r, 

IFd  .LT.4798)G0  TD  1 
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A. 


I i IW~  I |i 


O O O f3  O CJ  o o c 


I-I-l 

11  11=1 

97,  Wi^OT  ,7VTnn,r,Nij(  ID  , II 

97  -ORHAT(*  V99T  =*  , T 1 ?,  1 , » VTO®  =•  , FI  2 . 3 , • GN'J  =»,F12.3,*  II  *,I6) 

SUPPLY  MALF-<ir)HS  WMcsj  hj3T  ON  TAPE 

on  15  1=ILL,I1 
'IMr^^OL  (I) 

IF(NM.EO.l)  ro  TO  1*: 

[PlALPHAm.GT.l.M  GO  TO  13 
r-MM.EQ.2)  ALP1A(T>=0.07 
IF(NM,E0.3)  AL  P U ( Y)  =(J«  11 
IF(MM,E0.4)  ALP1A(T)=n.05 
IF<NM.t0.5)  A_P^A(T>-r0.1  j 
rF(YM,tO,6)  ALP^A^T)=^,(355 
t"(MM,t0.7»  ALPYA(Y)riO,o*5 

13  rF(ALPHA(I).LT.D.ri,OR, A.PHA(I).GT.l.O)  ALPHA < T) = 0 . 9 S 
15  CONTINUE 

calculate  the  ABSOpotion  COEFFICIENTS  AT  THE  'C  HOST  INTENSE  002  LINES 
ANT  AT  3 POINTS 


.OOF  BACK  PUNT  OON*T  NEED  TO  READ  TAPE 


301 

C 


312 
31  3 


314 

315 

C 

c 


309 


310 


OONTINUE 

NRITE<10,123)  i/l,/? 

»<RirE  (4, 120  '/l,\l? 

TATA  n£NS/1.0,l.  3,l.0'^-2,1.0E-3,  .5E-3,  .5E-2/ 

TO  512  JC=1,I1 

[F(GNU(  JC  I.GE.  >/l)  Gr<  tq  313 

:ONT INUE 

IHIN= JC 

TO  314  JC  = I'1IN,yi 
I^’(GNU(JC).GT.  v/2)  GO  TO  315 
TONTINUE 


IHAK= JC-1 


f W 


- AC=l.CE-20 

ILOOP=C 

JL00P=JL00P*1 

JPTS=0 

- A:  = , 3*F AC 

IHAX=IHAX^1 

GS'TOR=GNU(INAX) 

3ST0R=S<IMAX) 

1SiTOR  = MOL{I1AX) 

GMU( IhAX) *V? 

5 (I-AX)=l,C 
TO  37C  H=l,i 
iOL (INAX) =M 
TMEGA  {1,M»=  /1 
3HIN=FAC»0ENS(N) 

JCNT  = 1 

TO  32t  JC  = ININ,lriAX 
IF(NOL  (JC).ME.H»  GO  ■'O  320 
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ooo  ,v»oOf^, 


If fS< JC>.tT.Ml-4»  f,n  TO  3Sf 
DL»ST=OHEGA(JCNT,H) 

^ACl  = 2t.  ♦6N'J<  JC) 

IFFsFACl 
0tt£H  = IFF/20. 

IF(  <FAC1-IFF),GT.0.‘5>  OMi«*ONEM*0.0  5 

IO€LT  = 2e.  •<0W£W-0LA«:T» 

IFCIOELT-1)  320, 343, -^44 
T43  JCMT  = JCNT  + 1 
30  TO  319 

S44  3TE6A  (JCtlTn,’1>=Ot.fl''T^fl,  )5 
IF(IOElT-3)  345,31»i,317 
345  JC4r*JCt»T^2 
30  TO  319 

316  )»EGA  („iCNTf2,i)=OLa**T»o.lO 
JCMT= JCNT+3 

30  TO  310 

317  ONEGA ( JCNT»2,N) =.?• fOLAST»ONEW) 

ONEGA  (JCNT ♦ 3,1)  =0f‘"M-0.  0 5 
JCNT=JCNT+4 

319  ONEGA  (JCNT,1)=01Ef 

320  CONTINUE  ^ 

iCAtCm)=JC+»T  - ... 

320  Jfrs=JPTS+JCALC(i) 

5MUIIHA)0=G9T9R 
5 (IMAX)  = S5T0P 
10L(IMAX)=MST9R 
INAX=IMAX-1 

lF-(tilOOP.GE.60).ANn*  U®TS^»GT*240T»  GO  TO  371 
IFUJPTS.LT. 170  *ANn,  (FAC, GE. SLOWER)  ) GOTO  3 09 
IF< JPTS.LT.240)  GO  TO  321 
371  •AC=4.*FAC 
30  TO  389 


321  :01TINUE 
fCOUHT=0 
30  240  M=l,6 
HtAK=  JCALC(N) 

JN0NE=5 

If  I JtlAX.GT.  TMOME)  GO  TO  330 


WO  STRONG  LINES  IN  TlIS  BLOCK 

(MEGA  (4,M)«(0NEGA(9,N)f01EGA(3,H) )/2. 

3MEGA (2,M)=(01EGA«3,M>t09CGA(1,N>)/2, 

’RINT  339, JLOOP,FAC,S»ECIE(M) , (ONEGA (JJ,H) , JJ=1, J1AX) 
339  ^ORMAK//*  10  IlTENSf  LUES* , lOX  ,®JLOOR*»,  14  , 

IlOX,  •SHINa*,E10.3,l?X,*S»ECIESt*,  AiO/ 

1*  0MEGA«*,5F14.  T) 

30  TO  338 

STRONG  LINES  IN  THI‘^  9L03K 
330  GONTINUE 

»)MNT  S29,iNAX,  #LOn»,FAC,SPfCIE(N) 
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I 


325  -0?MftT  (//ItO,»  r.ftl  POIMTSt*,10X,*  JL03»  = *,  T3,5X, 

,t  *»S*1TN=*,P1C.  7,  nx,  61 ‘»/9X,*VflLUES  OF  OMEGA*) 

' >RINT  326,nMEGA(J),M)  , JJ  = 1,  JMAX) 

326  -0»MAT  (lOFl  .3) 

338  :OMTINUE 

.OSIC=.FALS-:, 

rF(M,e0,l,09.M.-’T.  LOGIC=.TRUE, 

i c 

C PRESSURE,  TEMOE^ATUO-^  LODP 

C 

c 

^0=1013, 00 
ra=296.oo 
')PT  = 0 

c 

C BRA'JCH  5ACK  POIJT  (“.T)  LlOP 

c 

350  -JPTrNPTH 
15  = 1 

:si= (TL-T  (NOT) ) / (Tr*T (NPT ) *0 . 6 946 ) 


i 

i 

i 


I 


fs 

G ROTATIONAL  OARTTTTO^'  FiJNGTIDN  IS  DEFINED  9EL0W 

C 

.IT  = 5QPT(T0/T(NPT)) 

:S'2(M)=T0/T(NPT) 

IF(loGIC)  CS2(M)  =0G'’(9)  *HT 
:fr=HT*P(NPT) /PC 

c 

C TEMPEPATURE  0"P“NO'‘Mr;F  0-  ALL  LINE  INTENSITIES  COlP'lTFD  HERE. 

C 

DO  23  I=ILL,I1 

IF (MOL  (I)  .NT.i)  GO  TO  2 ? 

50(r)=S(I)*G52(H)*FTO(-EDO(I)*CSl) 
aLPHAO(I)=ALPHA(I)*r6 
23  DOHTINUE 
C 

DO  555  JJ=l,jHAX 
/ = 0MEGA(  J 1) 
r 

:ari=c.o 

27  SUH1=C.0 

c 

C determine  indices  (T5  AND  16)  INDICATING  WHICH  SOECTR6L  LINES 

C ARE  TO  9£  USE')  IN  THC-  CA.CULATION  AT  FREQUENCY  V, 

C 

DO  '3  1=15,11 

IF  (V-e0UN0-GNU(I)»  '*9,29,33 
29  I5=r 

GO  TO  35 
33  CONTINUE 

15  = 11 
GO  TO  555 

39  DO  39  J=I5,T1 

1=  ( VTfiOUND-GNlH  J>>  TT.JT.SS 
37  T6=J-1 

GO  TO  43 
39  CONTINUE 
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I 


OOOC>  rt  O ! C^C^C90C»OrtOC^C^  ooo 


► 


:t)*»PUTC  THE  OPTICAL  n«-PTi  AND  TRANSilTTANCC  AT  F^iOUEHCY  V. 
3V»0.05 

43  30  45  Is  15, re 

IF  (HOL(I).^E 

AL«ALPHAO  (I) 

IF( ABS<7) «Gt« 

IF(AL.6T.0i95 

ARCTAN  FORmiA 

SUR1=S<HI)  /OV*  ( ATAM?  f Z+OY/2.  , At)-ATAN2  < Z-OT/2*,  ALM 
SO  TO  44 

.ORENTZ  FORMULA 

se  SUM1=S0(I)*AL* (.25/C (7-0Y/2.»**2^AL*AL) 

A^.50/  <Z*Z ♦AL^AL )♦.?*?/ nZ»-0V/2.  »**2+AL*AL)» 

44  :*Y1  = CAY14.S'JM1 

45  CWTINUE 

YSrOR(JJ)  =.318  3*CAYl*4(-i| 

555  SWTINVe 


, M)  GO  TO  45 


0*35)  TO  88 
) GO  TO  98 


- w i 


--i  A^i 


lQP/ 


30  557  JJ=1,JHAY 
557  5TD^ (NPT, JJ)=VSrOP(JJ) 

5*5  IF  (NPT.LT.NPTPTS)  GO  TO  350 


^RIITE  UROA  TA3LE 


1 


I 

I 


<«RITE(10,220)  S»ErTccH)  , JHAX 
<«»tTCtr,221)  SP€CIE<*t>,i'rAX  _ 

220  ~0?HATC1X,A4,I5) 
ttl  •0R5AT<A4,I5) 

DO  230  JJ=1,JMAX 
tCOUNT=ICOUNT^l 

IF  ( (ICOUNT.GE.?4C).ANO.(JHAX.NE.2))  GO  TO  230 
C 4RITEtl«,225)  OHfGA  ( J.l,  H>  f ISTORtNPT,  JJ»  , NPT=1  ,MPf  PT5) 

4RITE (4,226)  0MEGA(JJ,H) , (STOP (NPT, J J)  ,NPTsl,  NPTPTS) 

225  '0*MAT<F12.2,4<E12.6»/5(il2.6)  ) 

226  -0RHAT(F12.2,9(E12,6)  ) 

280  CO’Af IMtJE 

240  CONTINUE 

rt*»l»1f5L0CK 

IF(41.6E. Yrf2)  STOP  ?T 
VT  0 P= VI  ♦ YBLOCK  ♦ 501W  0 
IFIVTOP.GT.GNUdin  GO  T3  999 
V2SV14VBL0CK 
480T* Vl-BOUND 
- 30  TO  SO# 

C 

ERO 
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